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INTRODUCTION 


i . 


The  electron  be:im  technique  is  a well-established  method  for 
determining  species  concentrations  and  the  rotational  and  vibrational 
temperatures  of  molecular  nitrogen  in  low  density  flow  fields.  In  this 
method,  an  electron  beam  generator  produces  a narrow  column  of  electrons 
which  is  projected  across  the  flow  field  of  interest.  'Ihe  interaction 
of  the  accelerated  electrons  with  gas  particles  produces  a column  of 
radiation  which  is  nearly  coincident  with  the  beam  of  electrons.  Under 
aj  ropriate  conditions,  the  fluorescence  can  be  analyzed  spectroscopi- 
cally to  determine  the  concentration  and  temperatures  of  the  active 
species. 

The  predominant  radiation  due  to  excitation  of  air  is  the  first 
negative  system  of  the  ionized  nitrogen  molecule.  The  most  intense  band 
of  this  system  is  the  (0,0)  band  at  391^  A.  Much  experimental  work 
(summarized  in  Ref.  1)  has  been  done  with  the  first  negative  system  in 
nitrogen  and  air  to  verify  the  applicability  of  the  diagnostic  technique 
at  sufficiently  low  density  and  high  temperatures.  Details  of  the 
methods  can  be  found  in  Refs.  1-4. 

In  addition  to  radiation  in  the  No+  first  negative  system,  excita- 
tion of  air  flows  also  leads  to  observation  of  radiation  in  the  nitric 
oxide  y system.  The  > system  may  be  useful  for  vibrational  temperature 
and  number  density  measurements  at  wavelengths  below  approximately 
3044  A.  At  higher  wavelengths,  the  y bands  are  overlapped  by  bands  from 
II; > and  N>+  systems.  Additional  radiation  in  the  Feast  and  Ogawa  systems 
also  appear  at  wavelengths  between  5800  and  6100  A;  however,  these 
systems  are  unsuitable  for  electron  beam  diagnostics  because  of  their 
low  radiation  intensities  and  the  paucity  of  information  on  the 
exeitation-emisL .on  processes  which  lead  to  the  radiation.  (Details  of 
the  excitation-emission  analysis  for  the  NO  > system  are  given  in 
Ref,  5 and  typical  application  of  the  techniques  for  NO  are  given  in 
Ref.  6.) 

Electron  beam  excitation  of  air  flows  also  yields  radiation  in  the 
0->+  first  and  second  negative  systems.  The  first  negative  system  occurs 
at  wavelengths  near  6000  A and  is  not  suitable  for  electron  beam  measure- 
ments because  of  the  severe  overlap  of  vibrational  bands  within  the 
system.  'Ihe  0->+  second  negative  system  occurs  at  wavelengths  between 
2060  and  3230  A and  is  suitable  for  measurement  of  0;;  vibrational  tem- 
peratures. (The  basic  character  of  the  radiation  in  both  the  0 >+  first 
and  second  negative  systems  is  discussed  in  Ref.  7 and  applications  of 
the  measuring  techniques  are  given  in  Ref.  8.) 

To  apply  the  electron  beam  technique  for  measurement  of  rotational 
temperatures,  vibrational  temperatures , and  number  densities,  a detailed 
analysis  of  the  excitation-emission  process  which  leads  to  the  observed 
spectra  must  be  available.  It  is  coranon  to  assume  that  the  emission  is 
excited  by  collisions  between  beam  electrons  and  gas  molecules  in  their 
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ground  electronic  energy  state.  All  other  excitation  processes 
(excitation  by  secondary  electrons  and  low  energy  ions,  double  excita- 
tion, cascading,  etc.)  and  nil  selective  quenching  mechanisms  are 
ignored.  In  addition,  the  excitation  is  assumed  to  occur  so  rapidly 
that  no  perturbations  in  the  rotational  and  vibrational  population  dis- 
tributions in  the  ground  electronic  state  occur  during  the  excitation 
process;  hence,  the  Franck-Condon  principle  is  applied.  Rotational 
transition  probabilities  are  obtained  by  assuming  that  only  those 
transitions  occur  which  are  allowed  by  optical  selection  rules. 

At  temperatures  above  approximately  10*/  K.  and  densities  below  those 
corresponding  to  a pressure  of  1 torr  at  room  temperature,  the  assump- 
tion of  direct  excitation  and  application  of  optical  selection  rules 
lead  to  high  accuracy  in  density  and  temperature  measurements.  However, 
as  the  gas  density  is  increased,  both  collision  quenching  and  excitation 
by  secondary  particles  generated  by  the  primary  excitation  process  occur 
with  increasing  importance.  Both  of  these  mechanisms  invalidate  number 
density  measurements,  since  the  radiation  intensity  is  no  longer  a pre- 
dictable function  of  gas  density,  rotational  temperature,  and  vibra- 
tional temperature. 

There  is  much  interest  in  extending  the  density  range  over  which 
the  electron  beam  can  be  applied,  particularly  for  temperature  measure- 
ments. The  technique  is  attractive  for  a variety  of  flow  fields  because 
the  measurements  yield  temperatures  and,  in  some  cases,  actual  popula- 
tion distributions  without  causing  measurable  perturbations  in  the  flow 
field. 

Ashkenas0  reported  the  results  of  nitrogen  rotational  temperature 
measurements  made  in  a static  test  chamber  and  in  an  un-heated  free  jet 
with  equivalent  pressures  up  to  6.34  torr.  (Equivalent  pressure  is  de- 
fined as  the  static  pressure  obtained  frau  the  actual  gas  density  and 
a static  temperature  of  300  K.)  Ashkenas  found  that  the  error  in  measured 
rotational  temperature  increased  logarithmically  with  density.  Similar 
results  were  reported  by  Hunter.10  Maguire11  reported  rotational  tem- 
perature results  obtained  fretn  a free  jet  with  static  temperature  from 
00  to  300  K and  equivalent  pressures  up  to  10  torr.  These  free  jet  data 
show  that  excitation  by  low  energy  secondary  particles  is  not  a likely 
cause  for  errors  observed  in  rotational  temperature  measurements. 

Selective  quenching  of  the  rotational  energy  levels  was  believed  to  be 
the  cause  for  differences  between  measured  and  predicted  temperatures, 
if  such  selective  quenching  exists,  certain  rotational  energy  levels 
will  be  quenched  more  effectively  tnan  others,  so  that  the  excitation- 
emission  analysis  used  to  interpret  uhe  measured  intensities  in  terras 
of  rotational  and  vibrational  temperatures  cannct  be  applied. 

While  there  have  been  many  studies  of  the  range  of  applicability 
of  the  method  for  nitrogen  rotational  temperature  measurements  (see 
Muntz1),  little  attention  has  been  focused  on  similar  studies  or  the 
measurement  of  vibrational  temperatures.  The  techniques  have  been 
widely  applied  in  flows  with  high  vibrational  temperatures, >H»l  >_1H 


2 


but  the  gas  densities  have  been  in  the  range  where  the  techniques  are 
known  to  apply  with  acceptable  accuracy.  Harbour17  examined  vibrational 
band  intensity  ratios  in  low  temperature  N2  for  equivalent  pressures 
up  to  3 torr  and  found  large  variations  in  the  ratio  of  che  N2+(0,l)  and 
(l>2)  vibrational  bands.  However,  there  have  been  no  experiments  aimed 
at  evaluating  the  accuracy  with  which  vibrational  tempera  ^ure  measure- 
ments can  be  made  at  high  densities. 

The  present  studies  were  motivated  by  a requirement  to  examine  the 
accuracy  with  which  vibrational  temperatures  can  be  measured  with  an 
electron  beam  in  high  density  flow  fields  generated  by  facilities  with 
high  reservoir  temperatures . Calibration  experiments  were  conducted  in 
a static  test  chamber  and  in  an  arc-heated  wind  tunnel  facility.  The  equiv- 
alent pressures  ranged  from  0.12  torr  to  50  torr  with  static  temperatures 
from  300  K to  1000  K,  and  vibrational  temperatures  from  300  K to  3000  K. 
These  studies  were  conducted  to  assess  the  relative  inportunce  of  e: 'cita- 
tion by  low  energy  particles  and  excessive  collision  quenching  on  the 
accuracy  with  which  temperature  measurements  can  be  conducted  in  molecular 
nitrogen,  nitric  oxide,  and  molecular  oxygen. 


II.  THEORETICAL  ANALYSES 


A.  MOLECULAR  NITROGEN 


The  first  negative  system  in  molecular  nitrogen  is  due  to  radiative 
transitions  from  an  electronically  excited  state  of  N2+(N2+B?E)  to  the 
ion  ground  state  (N2+X?£) . The  excitation  results  from  collisions  of 
the  beam  electrons  with  nitrogen  molecules  in  their  ground  electronic 
state  (NaX1^).  With  the  assumption  of  direct  excitation  and  the  use  of 
optical  selection  rules,  the  excitation-emission  analysis  predicts  bands 
composed  only  of  a single  R -branch  and  a single  P-branch.1  The  P- branch 
lines  are  very  close  together  and  form  the  band  head,  while  the  R-branch 
lines  are  farther  apart  with  wavelength  separations  which  increase  with 
decreasing  wavelength.  The  intensity  of  a particular  line  can  be 
written  as4 


I - 


CSK-  v4q(v',v")  NN2  ^ q(Vo,v’) 

2K  ’ + 1 Qv(Tv)  £ QrTvo,  *r) 

[-  SijjalSs]  hK.KMr)/^  q(V.Y")  Kvii  (1) 

V ‘ / y" 


x exp 


where  Sg'  = K’  for  R-branch  lines  and  Sr'  = K'  + 1 for  P-branch  lines. 
The  vibrational  quantum  numbers  in  the  N2XJ£,  N2+B?E,  and  Na4^2!  states 
are  denoted  by  Vn,  v',  and  v",  respectively.  The  rotational  quantum 
number  in  the  N2  B?£  state  is  denoted  by  K'.  The  quantities  of  the 
form  q(v^,v^)  are  the  Franck-Con don  factors  for  the  -4 
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vibrational  transitions,  Qy  and  Qr  are  the  vibrational  and  rotational 
partition  functions,  and  v is  the  wave  number  of  the  line.  The  factor 
H(v0,K',Tr)  in  /riven  by1' 


(K ' + 1)  exp 


Bvn(K'  +•  l) (K  ' + 2 )hc' 
JcTr 


+ K exp 


BvnK ' (K ' - l)hc] 
kTR 


\2 ) 


The  rotator  constant,  Byo , is 

BVq  - 2.010  - .Ol8?(v0  + 1/2) 


(3) 


An  accurate  approximation  to  Eq.  1 is  obtained  by  neglecting  the 
dependency  of  the  rotator  constant  and  the  rotational  partition  func- 
tion on  the  degree  of  vibrational  excitation.  In  this  case,  Eq.  1 is 
written  as1 


where  0 Bhc/k  and  (G)  is  a known  function  of  Tp  and  K1  . The  rotational 
temperature  is  determined  by  measuring  the  relative  intensities  of  the 
rotational  lines  and  plotting  ^[i/Sk*  (G)v4]  versus  K'  (K’  + 1).  The  re- 
sult of  this  plot  is  a straight  line  with  a slope  equal  to  - 0/ Tr. 

Rotational  and  vibrational  temperatures  are  determined  with  a 
spectroscopic  instrument  yielding  line  and  band  intensities  in  the 
spectral  region  of  interest.  The  rotational  temperatures  are  usually 
obtained  from  the  rotational  lines  in  the  N2+(0,0)  band,  while  vibra- 
tional temperatures  are  determined  from  scans  cf  at  least  two  bands 
[typically  N?+(0,l)  and  N?+(l,2)].  For  vibrational  temperature  measure- 
ment, integration  of  the  band  intensities  should  be  used  to  obtain 
accurate  total  band  intensities. 

At  high  rotational  temperatures,  the  rotational  lines  corresponding 
to  high  rotational  quantum  numbers  overlap  the  low  rotational  quantum 
number  lines  of  the  next  band  in  the  sequence.  For  example,  rotational 
lines  in  the  (0,1)  band  with  rotational  quantum  numbers  greater  than 
24  overlap  the  P-bronch  of  the  (1,2)  band.  This  overlap  of  rotational 
lines  must  be  included  in  the  theoretical  analysis  to  allow  unambiguous 
interpretation  of  the  measured  band  intensities  in  terras  of  vibrational 
temperature . 

The  band  overlap  can  be  obtained  by  computing  the  detailed  intensity 
distribution  with  wavelength  and  will  depend  upon  the  resolution  of  the 
spectrometer  used.  This  requires  reasonable  knowledge  of  the  instrument 
function.  When  the  spectrometer  has  slit  widths  large  enough  that  the 
measured  line  width  is  much  greater  than  the  line  width  resulting  from 
broadening  processes,  the  instrument  function  can  be  obtained  by  assum- 
ing that  the  intensity  distribution  within  a line  is  flat-topped  as  the 
line  leaves  the  grating  of  the  spectrometer.  As  the  line  is  moved  past 
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the  exit  slit  (by  rotating  the  grating),  the  intensity  passing  through 
the  exit  slit  increases  linearly  with  wavelength  until  the  image  of  the 
entrance  slit  exactly  coincides  with  the  exit  slit.  For  equal  entrance 
and  exit  slit  widths,  further  grating  rotation  leads  to  a linear  decrease 
in  intensity  with  wavelength.  As  shown  in  Fig.  1,  the  ideal  instrument 
function  for  equal  entrance  and  exit  slit  widths  is  an  isosceles  tri- 
angle. The  ideal  line  profile  is  degraded  by  the  actual  intensity  dis- 
tribution within  the  image  of  the  entrance  slit,  astigmatism,  coma,  and 
improper  focusing.  Hie  effective  dispersion  is  obtained  with  the  ideal 
instrument  function  and  measurements  of  line  widths. 

Theoretical  estimates  of  band  profiles  within  the  electron  beam- 
induced  radiation  can  be  obtained  with  the  line  intensity  formulas  and 
the  effective  instrument  function  of  the  spectrometer.  A given  spectral 
region  is  divided  into  an  arbitrary  number  of  points  (typically  500)  at 
which  intensities  are  calculated.  At  each  of  these  points,  the  lines 
which  contribute  to  the  measured  intensity  are  determined  from  the 
wavelength  of  the  lines  and  the  line  width.  The  line  intensities  cal- 
culated from  Eq.  (l)  are  multiplied  by  the  instrument  function  (Fig.  l) 
and  summed  to  give  the  intensity  at  the  specific  point.  Hie  resulting 
intensity  distribution  is  integrated  to  give  the  total  band  intensities. 
The  calculations  are  easily  performed  with  a digital  computer. 

Typical  comparisons  of  measured  and  predicted  vibrational  band 
profiles  for  the  (0,1)  and  (1,2)  bands  of  the  N2+  first  negative  system 
at  vibrational  and  rotational  temperatures  of  300  K.  are  given  in  Fig.  2. 
Comparisons  of  the  measured  and  predicted  integrated  band  intensities 
are  given  in  Fig.  3-  The  only  significant  difference  between  the  experi- 
mental and  theoretical  band  profiles  occurs  at  the  band  head.  In  this 
region,  the  electro-optical  system  used  to  obtain  the  experimental  band 
profiles  did  not  respond  sufficiently  to  give  a true  measure  of  the 
intensity  distribution  within  this  high  wavelength  portion  of  the  (0,1) 
band  head.  Except  for  this  effect,  the  theoretical  model  gives  an 
excellent  representation  of  the  measured  band  profiles  for  the  Na+ 
first  negative  system. 

The  band  intensity  ratios  for  the  (0,1)  and  (1,2)  bands,  including 
the  effects  of  overlap  due  to  high  rotational  temperatures,  are  given 
in  Fig.  4.  It  is  to  be  noted  that  significant  errors  in  indicated  vi- 
brational temperature  occur  when  th-j  overlap  accompanying  high  rota- 
tional temperatures  is  ignored.  All  of  the  vibrational  temperature 
data  reported  here  were  obtained  from  Fig.  4 employing  the  measured 
rotational  temperatures  and  the  integrated  (0,1)  and  (1,2)  vibrational 
band  intensities. 


B.  NITRIC  OXIDE 

The  nitric  oxide  7 bands  result  from  a - ,vit  electronic 
transition  and  the  rotational  structure  is  quite  complicated.  In  the 
state,  the  electronic  orbital  angular  momentum  Is  zero  and  the 
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the  N2+  first  negative  (0,1)  and  ( l,2)0bar.ds; 
theoretical,  dispersion  = 14-7  A/mm; 

experimental,  PgQ  = 0.2  torr,  slit 

width  = 150  tim,  Tr  = Tv  = 300  K. 


INTEGRATED  INTENSITY 


Figure  3.  Theoretical  and  experimental  integrated  band  profiles 
for  the  N2+  first  negative  (0,1)  and  (J,2)  bands; 

theoretical,  dispersion  **  14.7  A/mm; 

- experimental,  PgQ  = 0.2  torr,  slit 
width  = 150  urn,  Tr  = Tv  * 300  K. 
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Figure  h.  Integrated  vibrational  band  intensity  ratios  for 
the  No+  first  negative  '^,x)  and  (1,2)  bands 


coupling  between  the  various  angular  momenta  can  be  assumed  to  follow 
ilund’s  case  (b).!“  The  rotational  energy  levels  in  the  ’upstate  con- 
sist of  two  components  corres ponding  to  J K ' l/2,  where  K denotes 
the  angular  momentum  vector,  apart  from  spin,  and  >f  is  the  total  angular 
momentum  vector. 

In  the  'n  ground  state,  the  electronic  orbital  angular  momentum  is 
non-zero  and  the  molecular  coupling  usually  can  be  assumed  to  follow 
either  Hund’s  case  (a)  or  (b).  In  case  (a),  the  rotational  energy 
levels  consist  of  two  components  denoted  by  iti/a  and  Vi/a>  where  the 
subscript  gives  the  values  of  the  quantum  number  associated  with  the 
sum  of  the  components  of  electronic  orbital  angular  momentum  and  spin 
in  the  direction  of  the  intemuclear  axis.  In  case  (b),  the  coupling 
is  the  same  as  for  >?  states  and  the  multiple  levels  correspond  to  J 
values  of  K I.  1/2.  It  is  often  found  that  the  molecular  coupling 
follows  case  (a)  at  lower  rotational  quantum  numbers  and  approaches 
case  (b)  at  higher  quantum  numbers.  While  such  is  the  case  for  the 
NO  y bands,  such  a transition  in  coupling  makes  calculation  of  wave- 
lengths difficult  at  intermediate  quantum  numbers.  Churchill19  gives 
term  values  which  apply  throughout  the  h’O  > bands,  circumventing  the 
need  to  choose  between  case  (a)  and  case  (b)  in  the  calculation  of  line 
positions . 

The  selection  rules  for  transitions  between  the  two  states  in  the 
NO  7 system  are  AT  - 0,  1 1.  Hence,  R,  P,  and  Q branches  exist  for 
transitions  from  each  multiplet  level  in  the  >'+  state  to  the  state. 
In  addition,  satellite  lines  appear  corresponding  to  transitions  for 
which  is  not  equal  to  AT. 

The  intensity  of  a rotational  line  in  the  NO  > system  resulting 
from  direct  excitation  of  nitric  oxide  with  an  electron  beam  is  given 
in  Ref.  5 as 


hcv‘tp(v'  ,v'')P(j  ' ,J")tp(v  1 ,J  'j 

v"  L j”  J 


(5) 


where 

rp(v  V')  c £ p(v>")|^  N(V",J")P(J',J'’)J  (6) 

where  P(J',J")  is  the  properly  normalized  Honl-London  factor  and  p(v',v") 
is  the  vibrational  transition  probability,  which  is  assumed  equal  to  a 
constant  times  the  Franck-Condon  factor.  The  Honl-London  factors  and 
Frnnck-Condon  factors  are  summarized  in  Ref.  5.  Ihe  term  values  used 
to  calculate  the  line  positions  are  given  in  Table  I. 
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Table  T.  Molecular  Constants  for  NO  v System 


«l/:2  «3/2 

Te 

0.0 

43965.7 

we 

1904.03  1903.68 

2371.1 

u'eXe 

13.97 

14.48 

WgYg 

-0.00120 

-0.28 

Be 

1.7046 

1.9930 

ae 

0.0178 

0.0198 

Be 

5 x 10■', 

6.0  x 10*" 

Be 

.... 

0.3  x lO*" 

Vibrational  Term  Values: 

G(v)  = we(v  + £-)  - uexe(v  + A)”  + weye(v  + ?)  ' 

Rotational  Term  Values; 

> O 

X1./2  > fl-3/2 

for  -«  3/ 
for  ■«!  / -> 

F(K)  - BVK(K  + 1)  - DyK^K  + 1)* 

Dy  De  + 0e(v  + +) 


F(.T)  - By 


l(J  + l):' . 

' J1’ 


1.0  X 0.5 


f 


Mj  + D9 

4j2 


124.2 

B„ 


- ')] 
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Ilic  numerical  methods  discussed  in  the  previous  section  were  used 
to  model  tl'.e  intensity  distribution  within  the  hands  of  the  NO  v system. 
Detailed  comparisons  of  the  calculated  wavelengths  of  the  lines  were 
made  with  the  data  given  in  Ref.  20.  The  present  theory  yields  line 
jHisitions  with  errors  no  greater  than  a few  wave  numbers.  A typical 
comparison  of  predicted  and  measured  vibrational  band  profiles  for  the 
(0,4)  and  (1,5)  bands  is  given  in  Kig.  l) . It  can  be  seen  that  the 
theoretical  model  slightly  underestimates  the  relative  intensities  of 
the  rotational  lines  in  the  tail  of  the  bands  and  that  tiie  predicted 
relative  magnitude  of  the  bands  within  a sequence  (av  constant)  are 
incorrect.  At,  room  temperature , the  theory  predicts  that  the  (l,v") 
bands  will  have  the  greatest  intensity  within  all  sequences.  This  re- 
sults from  a shift,  in  the  potential  curves  (energy  vs.  inter-nuclear 
separation)  between  the  NO  ground  state  and  the  r*  upper  state  of  the 
> system.  This  shift  causes  the  Kranck-Condon  factor  for  the 
v ’’  - 0-'  v’  r 1 transition  to  he  the  1 a rges I.  for  v"  * 0.  At  room  tempera- 
ture, populations  of  vibrational  levels  for  v"  '•  0 are  negligible  so 
that  the  excitation  transitions  primarily  populate  the  v'  1 level  in 
the  2>'+  state.  Hence,  the  theory  predicts  that  the  (l,v")  band  will 
dominate  each  sequence.  However,  the  experimental  measurements  clearly 
show  that  (0,v”)  bands  dominate  at  room  temperature. 

There  are  many  possible  causes  for  the  discrepancies  between  the 
theoretical  and  experimental  band  intensities . Mechanisms  such  us 
cascading,  selective  quenching,  self -absorption,  excitation  by  secondary 
electrons,  and  upper-state  relaxation  are  possible  causes.  Neither  cas- 
cading, selective  quenching,  self-absorption  nor  excitation  by  secondary 
electrons  are  likely  mechanisms  causing  inapplicability  of  the  direct 
excitation-emission  model.  Cascading  is  possible  because  of  the 
Kr'.',t  — A'  v+  transition  of  the  Keast  1 system,  which  populates  the  A">'+ 
level.  However,  the  spectral  studies  of  Ref.  5 show  that  the  radiation 
in  the  Keast  system  is  extremely  weak,  indicating  that  relatively  few 
molecules  are  excited  to  the  bl  "+  state. 

Selective  quenching  of  the  vibrational  energy  levels  within  the 
: 5.  state  is  also  not  a likely  cause  for  the  discrepancies  between  the 
theoretical  and  experimental  relative  band  intensities.  As  shown  in 
Ref.  5,  the  NO  > system  has  a relatively  high  self-quenching  pressure 
of  2.3?  torr  and  ratios  of  band  intensities  are  constant  over  a range 
in  pressure  at  least  from  0.2  to  3.0  torr.  In  addition,  the  discrep- 
ancies are  observed  at  low  pressures  where  quenching  can  be  ignored. 

As  discussed  in  Ref.  5,  the  NO  > system  is  a resonant  transition 
so  that  radiation  emitted  in  the  region  of  the  electron  beam  may  be 
absorbed  by  NO  molecules  in  the  surrounding  gas.  This  re-ubsorption 
can  lead  to  induced  radiative  transit 'o ns  and  volume  fluorescence 
throughout  the  gas.  However,  the  total  light  intensity  due  to  the 
volume  fluorescence  and  induced  emission  can  be  only  as  large  as  the 
absorbed  light  intensity.  Hurther,  the  fluorescence  and  induced  emis- 
sion are  radiated  spherically,  sc  that  only  a small  percentage  of  their 
total  radiation  can  enter  the  optical  system.  Heine,  fluorescence  and 
induced  emission  can  be  ignored. 
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(1,5) 


Figure  5.  Theoretical  and  experimental  band  profiles 
for  the  NO  y (0,4)  and  (1,5)  bands; 

••••’•  tneoretical,  direcj  excitation, 
dispersion  * 14. 7 A/ram 

experimental,  Pj.;q  r-  0.3  torr,  slit 

width  = 150  urn,  Tr  = Tv  = 300  K. 
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Detailed  calculations  of  the  absorption  to  be  expected  in  the  NO  y 
system  are  fivers  in  Ref.  f;  where  it  is  shown  that  absorption  can  be 
neglected  for  NO  number  density-path  lengths  products  less  than 
10'“  cm"  when  the  emission  transition  is. to  vibrational  energy  levels 
other  than  the  ground  state  (v"  0).  In  the  present  studies,  the 

number  density-path  length  product  was  always  less  than  10^  cm-'’ , so 
that  the  effects  of  absorption  can  be  neglected. 


Excitation  by  low  energy  secondary  electrons  also  is  not  likely. 

•'hen  the  excitation  is  dominated  by  secondary  electrons,  the  intensity 
of  the  radiation  varies  quadratically  with  density  at  densities  low 
enough  to  avoid  significant  quenching. However,  the  results  of  Ref.  5 
show  that  the  variations  of  the  band  intensities  with  density  are  linear 
at  equivalent  pressure  below  0.20  torr.  Hence,  excitation  by  low  energy 
secondary  electrons  cannot  be  a predominant  mechanism. 

The  apparent  inapplicability  of  the  direct  excitation  model  for 
the  NO  y system  is  most  likely  related  to  the  polar  nature  of  the 
molecule.  Since  vibration-rotation  transitions  can  occur  optically 
within  an  electronic  energy  state,  redistribution  of  the  energy  can 
occur  both  in  the  X rt  state  in  the  early  stages  of  excitation  by  the 
high  speed  electrons  and  in  the  Ap£+  state  before  the  electronic  de- 
excitation transition  has  a chance  to  take  place.  Since  the  Born- 
Oppenheimer  approximation  is  expected  to  be  valid  for  excitation  by 
beam  electrons,  redistribution  of  the  energy  within  the  Xp«  state  prior 
to  the  ’x  - ~T.  electronic  transition  is  not  likely.  However,  redis- 
tribution of  the  energy  in  the  ?E+  state  due  to  vibration-rotation 
transitions  appears  quite  likely  and  may  influence  the  relative  band 
Intensities.  In  this  case,  the  vibration-rotation  transitions  must 
occur  within  the  lifetime  of  the  PE+  state. 

The  influence  of  vibration-rotation  transitions  within  the  p£+ 
state  prior  to  the  PE+  — > electronic  transition  can  be  estimated  by 

Including  upper-state  vibration-rotation  transitions  in  the  excitation- 
emission  model.  This  is  accomplished  by  assuming  that  the  NO  molecule 
in  the  I;  state  can  be  represented  as  a harmonic  oscillator.  In  this 
case,  vibrational  transitions  can  occur  only  between  adjacent  vibrational 
energy  levels  and  the  transition  probabilities  are  proportional  to  the 
vibrational  quantum  number  of  the  initial  state  in  the  transition.5,1 
Hence,  the  vibrational  transition  probabilities  can  be  written  as 

p*(v',v'  i 1)  = k10v'  (7) 

where  klo  is  the  probability  for  the  1 -*  0 vibrational  transition.  The 
allowed  rotational  transitions  are  those  for  which  AK  = ± 1 and  the 
correspanding  transition  probabilities  are  given  by  the  usual  Hbnl- 
:.ondon  factors. 


The  transitions  in  the  upper-state  are  included  in  the  excitation- 
emission  model  by  assuming  that  only  de-excitation  transitions  occur  in 
Hie  ’a+  state.  Hence,  all  collisional  effects  and  absorption  within 
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the  state  are  ignored.  With  these  assumptions,  the  rate  of  popula- 
tion of  a particular  (v',K')  level  in  the  state  can  be  given  by 

<p(v  1 »K’)  --  £££  N(vV")P(J‘,J")]  P(v',v") 

+ Ck(v‘  + l)k10[N(v'  + 1,K ' + 1)PRVR  + N(v'  + 1»K’  " 1)ppvp3 

(8) 


where  PR  and  Pp  are  the  Hbnl-London  factors  for  the  R-branch  and  P-branch 
rotational  transitions  and  C-p  and  Ck  are  constants  independent  of  the 
quantum  numbers.  The  first  term  on  the  right  hand  side  of  Eq.  (8) 
accounts  for  the  usual  direct  excitation  from  the  state  while  the 
second  term  accounts  for  vibrational  de-excitation  transitions  within 
the  state.  Note  that  only  P and  R branches  exist  within  the  ,?Sf 
state  and  that  the  splitting  accompanying  the  doublet  structure  is 
ignored.  The  rotational  and  vibrational  quantum  numbers  in  the  X’n 
state  are  donated  by  J"  and  v",  respectively,  and  the  sum  over  J”  in- 
cludes all  branches  in  the  composite  model  discussed  above. 


The  total  rate  of  depopulation  of  the  (v  ’ ,K ' ) level  is  assumed  to 
result  only  from  optical  transitions  and  can  be  written  as 


*(v',K')  = N(v',K‘)  fca£;  Ep(V,v-)P(jVVv.iT..J.1j. 

L v J”  * 

+ eKv’kl0(PRv^  + Ppvg)J 


(9) 


In  the  steady  state,  the  total  rate  of  excitation  must  be  equal  to  the 
rate  of  de-excitation  so  that  <p(v',K’)  - ij/(v',K').  Hence,  frcm 
Eqs.  (8)  and  (9)>  the  population  of  the  N(v',K')  states  are  given  by 

N(v',K')  - Jc^  £|"£  N(V',J")P(J'  ,J")J  p(v'?v”) 

+ CR(v'  + l)klu[N(v’  + 1,K'  + l)PRvg  t N(v’  + 1,K*  - l)PpVpl  J/' 

Ca  L E^v'>v">P(J'»J'')yv’,v",J’,T"  + Ckv’k10lPRv-  + Ppv|] 

1 V"J”  (10)1 


Examination  of  Eq.  (10)  shows  that  it  represents  a hierarchy  of 
equations  for  the  populations  of  all  the  vibrational  levels  in  the 
state.  To  obtain  e solution,  n maximum  value  for  the  vibrational  quantum 
number  in  the  r?T?  must  be  chosen.  Beezi^°  shows  that  the  NO  molecule 
dissociates  into  N(‘'S)  and  0(;,P)  at  an  energy  slightly  higher  than  that 
for  the  v'  = 3 level.  In  addition,  no  vibrational  bands  for  v'  > 3 
have  been  reported  for  the  > system.5 nonce,  it  is  reasonable  to 
terminate  the  vibrational  quantum  number  at  3 in  Eq.  (10).  For  v*  » 3 
there  is  no  population  due  to  transitions  from  higher  vibrational  levels 
so  that  Eq.  (10)  can  be  written  as 


rr-  N(3.K'  ) 

°<P 


T N(vV")P(J  ',J")p(v*,v") 


£1  p(v',v")P(J' ,J")v3  + 3(Ck/Ca)ku,[PRvg  + Ppvg] 
V"  J" 


(11) 


Assuming  a Boltzmann  distribution  for  the  X3n  state  allows  relative 
values  for  N(v",J")  to  be  obtained  so  that  Eq.  (ll)  can  be  solved  for 
the  relative  value  of  N(3>K’).  With  this  value,  the  population  of  the 
remaining  vibrational  energy  levels  can  be  solved  sequentially  for  de- 
creasing values  of  v’ . Once  the  populations  are  obtained,  the  relative 
intensities  of  the  rotational  lines  in  the  y system  can  be  obtained  from 

I - av4p(v',v'')P(J‘ ,J”)N(v',K' ) (12) 

where  a is  a constant  of  proportionality. 

Explicit  solution  of  Eq.  (11 ) requires  that  the  relative  value  of 
the  vibrational  transition  probability  (Ck/Ca)kl0  be  available.  The 
term  Ck  relates  k10  to  the  absolute  vibrational  transitional  probability, 
which  can  be  obtained  with  good  accuracy  for  the  assumed  harmonic  oscil- 
lator model.  However,  Ca  includes,  in  part,  the  electronic  transition 
moment;  its  evaluation  requires  knowledge  of  the  absolute  electronic 
transition  probabilities  for  the  bands  in  the  y system. 

Since  little  precision  in  the  absolute  electronic  transition 
probabilities  is  available,  an  effective  value  for  (Ck/Ca)kXo  was  deter- 
mined by  comparing  theoretical  and  experimental  band  profiles  at  roan 
temperature.  The  value  of (Ck/Ca )kl0  was  adjusted  until  the  integrated 
relative  intensities  nf  the  (0,4)  and  (1,5)  bands  obtained  from  the  theory 
agreed  with  those  measured  at  pressures  below  0.3  torr  at  room  temperature. 
The  appropriate  value  for  (ck/cn)ki0  with  the  non-dimensionalizing  scheme 
used  in  Eqs.  (8-10)  is  0.08. 

Equations  (10-12)  were  used  to  compute  the  variation  of  the  ratio 
of  the  intensities  of  the  (0,4)  and  (1,5)  bands  with  rotational  and 
vibrational  temperatures.  In  addition,  the  band  intensity  ratio  obtained 
with  simple  direct  excitation  (klc)  - 0)  was  computed.  The  results  are 
given  in  Fig.  6.  It  is  noteworthy  that  for  the  spectrometer  resolution 
used  in  these  calculations  there  is  no  significant  overlap  of  the  (0,4) 
and  (1,5)  bands  for  rotational  temperatures  from  300  K to  1000  K. 

Hence,  tne  band  intensity  ratios  for  the  direct  excitation  model  are 
independent  of  rotational  temperature  over  this  temperature  range. 

However,  when  vibration-rotation  transitions  are  allowed  to  occur  in  the 
state,  the  band  intensity  ratio  depends  upon  both  the  rotational  and 
vibrational  temperatures.  For  rotational  temperatures  up  to  1000  K,  the 
theory  results  in  negligiDle  overlap  of  the  high  rotational  quantum 
number  lines  of  one  band  with  the  head  of  the  next  band  in  the  sequence. 
Hence,  the  variation  of  band  intensity  ratio  with  rotational  temperature 
implies  that  the  transitions  within  the  state  slter  both  the  rota- 
tional ana  vibrational  population  distributions  prior  to  the  emission 
uransitions . 
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In  practical  applications  of  an  electron  beam  to  measuring 
vibrational  temperature,  the  relative  sensitivities  of  the  optical 
system  are  obtained  by  measuring  the  relative  band  intensities  at  a 
known  condition  (usually  room  temperature).  The  resulting  scale  factor 
is  used  to  convert  the  measured  band  intensities  obtained  in  the  unknown 
environment  to  vibrational  temperature.  In  this  manner,  the  technique 
always  gives  the  correct  temperature  at  the  known  condition.  Hie  effect 
is  to  non-dlmensionalize  the  curves  of  Fig.  6 by  their  values  at  the 
known  condition.  The  band  intensity  ratio  curves  of  Fig.  6 non- 
dimensionalized  by  their  respective  values  at  rotational  and  vibrational 
temperatures  of  300  K are  given  in  Fig.  7*  Ibis  displays  the  relative 
importance  of  vibrational  relaxation  in  the  ?P+  state  and  it  is  seen 
that  the  non-dimensional  variations  of  the  band  intensity  ratios  with 
vibrational  temperature  are  altered  only  slightly  by  the  vibrational 
relaxation.  Hence,  the  main  effect  of  the  relaxation  is  to  change  the 
relative  values  of  tne  band  intensities  and  to  introduce  a rotational 
temperature  dependence . There  is  little  overall  effect  on  the  sensi- 
tivity of  band  intensity  ratios  to  vibrational  temperature  when  pre-run 
scaling  of  relative  band  intensities  is  employed. 

Comparing  Figs.  4 and  6,  it  can  be  seen  that  the  band  intensity 
ratios  for  the  NO  y system  are  much  less  sensitive  to  vibrational  tem- 
perature than  are  the  ratios  for  the  Nr>+  first  negative  system.  Hence, 
extreme  accuracy  is  required  to  obtain  reasonable  values  for  the  nitric 
oxide  vibrational  temperature.  An  error  in  measured  intensity  ratio  as 
small  as  5$  results  in  a vibrational  temperature  error  near  20%. 

The  applicability  of  the  excitation-emission  model  employing 
vibrational  relaxation  in  the  state  can  be  determined  only  by  de- 
tailed comparisons  of  theoretical  and  experimental  results  at  elevated 
vibrational  temperature.  There  is  no  guarantee  that  the  model  is  cor- 
rect since  a one-point  comparison  was  used  to  establish  a single  unknown 
relative  transition  probability.  Certain  comparisons  of  the  results 
from  this  model  with  band  profiles  measured  at  moderate  vibrational 
temperatures  are  given  in  the  next  section  of  this  report. 


C.  QUENCHING  MECHANISMS 

The  analyses  of  the  preceding  sections  apply  only  when  collision 
quenching  can  be  ignored.  A more  general  analysis  can  be  formula oed  to 
include  other  excitation/de-excitation  mechanisms.  To  allow  specific 
equations  to  be  formulated,  the  following  analysis  is  presented  for 
molecular  nitrogen  in  an  arbitrary  gas  mixture.  However,  the  results 
are  applicable  in  principle  to  any  radiating  species. 

The  important  excitation  processes  are  assumed  to  include  only 
excitation  by  primary  (beam)  electrons  and  secondary  electrons  liberated 
by  the  bombardment  of  ground  state  particles  by  primary  electrons. 

These  processes  are  represented  (for  N2)  by: 


VIBRATIONAL  TEMPERATURE  (K) 


(0,4)/(l  ,5)  Bond  Intensity  Rotio 


Figure  7.  Nondimensional  band  intensity  ratios  for  the  NO  7 (0,4) 
and  (1,5)  bands.  Ratios  nondimensionalized  by  their 
values  at  Tp  Tv  300  K 
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(13) 


N->(X  1.v0  + ep  Mvgl>  N :+(B;-v)  + ep  + es 

and 

Np(Xjz:)  + es  Jill— iL  N2+(B:’e)  + 2es  (l4) 

where  ep  and  es  denote,  respectively,  primary  and  secondary  electrons 
and  Qo(ve)  and  Qs(vs)  are  the  cross-sections  for  excitation  by  primary 
electrons  with  velocity  ve  and  secondary  electrons  with  velocity  vs, 
respectively. 

Secondary  electrons  are  ejected  from  the  parent  molecules  in  random 
directions  and  with  a range  of  energy.  Their  distribution  must  be  cal- 
culated by  balancing  the  effects  of  diffusion,  collisions,  and  space 
charge.  An  analysis  sufficient  to  describe  the  local  secondary  electron 
concentration  has  been  given  previously7  and  yields  a rate  of  excitation 
given  by 

nsvsng  £ xjQjT(ve) 

i 

XjQjs(Vs)  + xs/ng 

j 

where  ns  and  rig  are  the  secondary  electron  and  total  gas  number 
densities  and  Xp  is  the  number  density  fraction  cf  N2  molecules  in  the 
mixture.  The  total  cross-section  for  excitation  of  species  j by 
secondary  electrons  is  given  by  QjS(vs)»  where  Xj  is  the  number  density 
fraction  of  species  j,  and  Xg  is  related  to  the  diffusion  coefficient. 
The  term  Qj‘p(ve)  is  an  effective  ionization  cross-section  for  production 
of  secondary  electrons  in  the  energy  range  necessary  to  excite  N^X1?,) 
to  the  Na^B^x)  state. 

De-excitation  of  the  N2+{ByX)  is  assumed  to  result  from  radiative 
transitions  and  collision  quenching.  These  are  represented  by 

n3+(b :>r.)  . A— » n^Cx^)  + hv  (16) 

0Ild  Q, 

n2+(b:j5;)  + Nj  — ^ n2+(x"’>:.)  + Nj  (17) 

where  Nj  denotes  the  concentration  of  species  j,  Qjc  is  an  effective 
quenching  cross-section,  and  A is  the  sum  of  all  the  optical  transition 
probabilities. 

In  the  steady  state,  Eqs.  (13-17)  yield  the  population  of  the 
Na^B^X.)  state  as7 
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N?+(B'”>')  ] 


nevengQo(ve)XNo 
A[1  + V*  Z XG^c] 

Qs(vs)XNanengve  £ XjQjT(ve) 

+ ng/A  ^ ^j^js^vs)  + Xg/ngJ 


(18) 


Fo:*  the  present  discussion,  we  ignore  the  detailed  spectral 
structure  of  the  radiation  and  write  the  radiation  intensity  as 

I --  hcv'1A[Ng+(B£>S)] 

combining  Eqs.  (l8)  and  (19)  gives 


I 


angJQ0(ve)xMg  + angJ^S(vS)XNo 

i+L(nj/ny)  i+L(nj/nj‘) 

«j  J 


H xjQjT(ve) 

_J 

£ XjQjS(vs)  + As 


(19) 


(20) 


where  J is  the  beam  current  and  a is  a constant  depending,  in  part,  on 
geometrical  factors.  The  quantity  nj ’ is  the  quenching  density  of 
species  j and  in  a single  species  gas  gives  the  density  at  which  the 
intensity  is  reduced  to  one-half  the  value  obtained  in  the  absence  of 
collision  quenching. 


If  excitation  by  secondary  electrons  can  be  ignored,  the  intensity 
expression  degenerates  to  a Stern-Volmer'''3  equation  of  the  form 


I - 


aingXN£ 


+ £ (W* 


(21) 


At  high  densities.-  the  terms  in  the  summation  of  Eq.  (21)  dominate 
the  denominator  and  the  intensity  can  be  written  as 


I = 


(22) 


Note  that  if  one  term  of  the  summation  dominates,  Eq.  (22)  becomes 

I = ainj'XN2/Xj  (23) 

while  for  a single  species  gas,  the  intensity  is 

I - aj,n'  (24) 
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where  n’  is  the  self- quenching  density.  Under  the  conditions  lending 
to  either  Eq.  (23)  or  (2L),  the  radiation  intensity  is  independent  of 
the  absolute  number  density  and  becomes  n function  of  temperature 
(through  the  unique  dependence  of  n:'  on  temperature)  and  the  relative 
concentrations  of  N;>  and  the  quenching  species. 

It  is  clear  that  in  the  presence  of  severe  collision  quenching, 
measurements  of  the  absolute  number  density  cannot  be  made  with  an 
electron  beam.  It  may  be  possible  to  determine  relative  concentrations, 
but  only  if  the  detailed  dependencies  of  the  quenching  number  densities 
( n j ' ) on  temperature  are  known. 

To  evaluate  the  effect  of  selective  quenching  on  temperature 
measurements,  Eq.  (21)  is  written  as 


I - 


1 + 


—la 

L (nj/n,j' ) 

<3 


(25) 


where  I0  is  the  intensity  which  would  be  observed  if  there  were  no 
quenching  and  is  given  by  Eq.  (l).  Rotational  and  vibrational  tempera- 
tures are  usually  obtained  from  measurements  of  the  relative  intensities 
of  rotational  lines  (or  groups  of  lines).  If  the  quenching  is  not 
selective,  the  denominator  of  Eq.  (25)  will  have  the  same  value  for  all 
lines  and  will  not  alter  an  intensity  ratio.  In  other  words,  rotational 
and  vibrational  temperatures  can  be  measured  in  the  presence  of  quench- 
ing if  the  ny  values  do  not  depend  upon  the  rotational  and  vibrational 
quantum  numbers  in  the  N;>+(B?>:.)  state.  The  experimental  studies  de- 
scribed in  the  following  sections  of  this  report  were  directed,  in  part, 
at  examining  the  accuracy  with  which  vibrational  temperatures  can  be 
measured  when  the  intensities  resulting  from  electron  beam  excitation 
must  be  represented  in  the  form  of  Eq.  (25). 


III.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

A.  TEST  FACILITIES 

1.  Static  Chamber 


A series  of  experiments  was  conducted  in  a static  test  chamber 
shown  schematically  in  Fig.  8.  ‘Ihe  chamber  consists  of  a cadmium 
plated  steel  cylinder  18  inches  in  diameter  and  18  inches  long.  Vacuum 
conditions  in  che  chamber  are  maintained  by  a mechanical  vacuum  pump. 
Pressures  are  measured  with  a thermocouple  gage,  Wallace  and  Tiernan 
absolute  pressure  gages,  and  variable  reluctance  pressure  transducers. 
The  electron  beam  is  projected  across  the  tank  in  a direction  perpendi- 
cular to  the  tank  axis.  The  drift  tube  of  the  electron  beam  generator 
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Figure  8.  Static  test  chamber  schematic 
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is  electrically  isolated  fran  the  chamber  so  that  the  entire  chamber 
serves  as  the  beam  collecting  device.  A specially  designed  electron 
trap  is  used  to  reduce  reflection  of  beam  electrons  back  into  the  test 
chamber . 

The  electron  beam  generator  employed  for  these  tests  utilized  a 
duoplasmatron  electron  so’irce  and  is  described  in  detail  in  Ref.  5* 

The  nominal  beam  voltage  was  20  kV  and  the  beam  current  was  approximately 
2 mA.  Electron  beam  currents  were  measured  with  a standard  miliiammeter 
and  were  recorded  by  measuring  the  voltage  drop  across  the  meter  resis- 
tor. 


Experiments  were  conducted  with  nitrogen  and  air  at  static  pressures 
up  to  50  torr.  Pressures  were  maintained  by  balancing  the  pumping  rate 
of  the  vacuum  system  with  gas  injection  controlled  by  a precision  needle 
valve.  For  most  experiments,  the  gas  was  injected  at  the  top  of  the 
tank  and  the  vacuum  pump  was  connected  near  the  bottom  of  the  tank  in  an 
attempt  to  force  direct  flow  across  the  region  of  observation. 

At  static  pressures  greater  than  approximately  10  torr,  severe 
beam-plasma  interactions  were  observed.  These  were  manifested  as  large 
fluctuations  in  measured  beam  current.  In  most  cases,  negative  beam 
currents  were  observed.  The  erroneous  beam  current  signals  resulted 
from  potential  differences  between  the  insulated  drift  tube  and  the  test 
chamber,  so  that  the  drift  tube  collected  charged  particles.  This  situ- 
ation was  alleviated  by  connecting  a vacuum  tube  voltmeter  and  a dc 
power  supply  between  the  drift  tube  and  the  static  chamber.  The  power 
supply  potential  was  adjusted  to  maximize  the  beam  current  indication. 
However,  the  precise  value  of  the  beam  current  under  these  conditions 
was  unknown. 

2.  Arc -Heated  Wind  Tunnel 


An  arc  heated  wind  tunnel  facility  was  employed  to  provide  air 
flows  at  known  thermodynamic  conditions.  The  facility  is  snown  sche- 
matically in  Fig.  9 and  consists  of  a direct  current  arc  heater,  a 
settling  chamber,  a conical  convergent-divergent  nozzle,  a free  jet  test 
cabin,  and  a pressure  recovery  system.  The  settling  chamber  is  employed 
to  reduce  spatial  and  temporal  fluctuations  in  gas  properties.  In 
addition,  several  bare-tube  heat  exchangers  consisting  of  water-cooled 
copper  tubing  were  inserted  in  the  flow  to  reduce  the  flow  stagnation 
temperature.  In  certain  cases,  unheated  air  was  admixed  with  arc  heated 
air  in  the  settling  chamber  upstream  of  the  heat  exchanger  section. 

With  this  arrangement,  the  reservoir  temperature  could  be  varied  from 
approximately  1400  K to  WOO  K at  a nominal  reservoir  pressure  of 
1.0  atmosphere. 

Two  conical  nozzles  were  employed.  The  first  nozzle  (Nl)  had  an 
expansion  half  angle  of  3.5°  and  delivered  a flow  with  a nominal  Mach 
number  of  2.0  and  a static  pressure  of  60  torr.  The  vacuum  pumping 
station  maintained  a pressure  in  the  free  jet  test  cabin  of  approximately 
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Figure  9«  Arc-tunnel  schematic 


1 torr;  hence,  u "barrel-shock"  system  was  generated  by  the  highly  under- 
expanded  nozzle  flow  field.  The  electron  boom  was  projected  aeronn  the 
flow  field  at  various  distances  downs tre urn  of  the  nozzle  exit  to  provide 
test  regions  with  static  pressures  varying  from  1.5  to  60  torr  and  flow 
Mach  numbers  from  2.0  to  4.25. 

The  second  nozzle  (N2)  had  an  exit  expansion  angle  of  15°  und  de- 
livered a flow  with  u nominal  Mach  number  of  7.0.  The  nozzle  exit  dia- 
meter was  4.0  inches.  Since  the  nozzle  exit  and  test  cabin  presoures 
were  nearly  matched  for  this  nozzle,  the  axial  and  radial  gradients  in 
flow  properties  were  much  less  than  those  for  the  first,  nozzle.  The 
electron  beam  was  projected  across  this  flow  1.0  inch  downstream  of 
the  nozzle  exit. 

Wind  tunnel  reservoir  conditions  were  determined  by  measuring  the 
reservoir  pressure  and  the  total  mass  flow  rate  and  assuming  complete 
thermodynamic  and  chemical  equilibrium  from  the  reservoir  to  the  nozzle 
throat.  Under  these  conditions,  the  mass  flow  rate  and  reservoir  pres- 
sure are  uniquely  related  to  the  reservoir  enthalpy. 

In  nozzle  Nl,  the  electron  beam  was  generated  by  a commercially 
available  19AYP4  television  electron  gun.  Beam  currents  of  approximately 
0.5  mA  at  voltages  near  15  kV  were  employed.  A detailed  description  of 
this  beam  generator  and  some  of  the  special  operating  procedures  required 
are  discussed  in  Ref.  24.  For  nozzle  N2,  the  electron  beam  was  produced 
by  a duoplasma tron  generator  described  in  Ref.  5.  This  generator  pro- 
vided beams  with  currents  up  to  10  mA  at  voltages  near  20  kV.  For  both 
beam  generators,  the  entire  generator  system  was  electrically  isolated 
from  the  wind  tunnel,  so  that  no  special  beam  collecting  devices  were 
required. 


B.  INSTRUMENTATION 
1.  General 


The  wind  tunnel  was  instrumented  to  record  static  and  total 
pressures,  mass  flows,  and  inlet  gas  temperatures.  Reservoir  pressures, 
test  section  pressures,  and  pitot  pressures  were  obtained  with  variable 
reluctance  differential  pressure  transducers.  All  pressure  transducers 
were  connected  to  n vacuum  reference  which  was  held  at  approximately 
.005  torr.  The  transducers  were  calibrated  against  a silicone  oil 
micromanometer . 

Primary  (arc-heated)  and  secondary  (cold  gas  injection)  muss  flow 
rutea  were  monitored  by  measuring  the  differential  pressure  drop  across 
calibrated  orifice  pistes.  Ihe  pressure  drops  were  measured  with 
variable  reluctance  pressure  transducers . 


2.  Optical  System 


Rotational  and  vibrational  spectra  were  obtained  with  a Jarrell- 
Ash  O.^  Meter  Ebert  scanning  spectrometer.  'L'he  grating  was  ruled  with 
30,000  lines/inch  yielding  a dispersion  of  16  A per  millimeter  and  a 
maximum  resolution  of  0.2  A.  An  uncooled  EMI  6256  photomultiplier  was 
used  with  the  spectrometer.  A schematic  of  the  optical  system  is  given 
in  Fig . 10 . 

Hie  image  of  the  electron  beam  was  focused  on  the  entrance  slit  of 
the  spectrometer  by  a 13  in.  focal  length,  2-1/2  in.  diameter  quartz 
lens.  The  electron  beam  was  oriented  perpendicular  to  the  entrance 
slit.  Photocurrents  were  measured  and  amplified  by  a Keithley  Model  4l? 
picoammeter.  l'he  output  of  the  amplifier  was  entered  into  the  signal 
conditioning  system  of  the  AARL  Digital  Computer  and  Data  Acquisition 
System. 


The  electron  beam  current  was  continually  monitored  by  a second 
photomultiplier-spectrometer  system.  Radiation  from  the  electron  beam 
was  focused  on  the  entrance  slit  of  a 0.25  meter  Jarrell-Ash  spectrometer 
set  tc  view  the  head  of  (0,0)  band  of  the  N->+  first  negative  system. 
During  all  spectral  scans,  the  output  from  the  0.5  meter  spectrometer 
was  divided  by  that  from  the  0.25  meter  spectrometer  to  account  for 
variations  in  beam  current  and  gas  density  during  the  wavelength  scan. 

All  data  reported  here  were  obtained  in  this  manner. 

3.  Data  Acquisition  System 

The  experimental  data  were  collected  with  the  AARL  Digital 
Computer  and  Data  Acquisition  System.  This  system  is  employed  for  on- 
line data  acquisition,  facility  control  and  data  reduction,  as  well  as 
for  conventional  digital  numerical  processing.  For  these  tests,  a 
medium  speed,  auto-ranging  amplifier  and  analog- to-digital  (A/d)  con- 
verter system  was  employed.  The  A/D  converter  resolves  each  data  signal 
to  3.1  digital  bits  plus  a sign  bit.  Hence,  the  accuracy  of  the  con- 
verter is  + .0489%  of  full  3cale  on  all  gain  ranges.  Full  scale  voltage 
is  adjustable  from  .005  to  1.0  V in  8 steps.  The  appropriate  gain  is 
automatically  selected  by  the  computer  program  operating  the  system  to 
maximize  data  accuracy. 

Data  were  obtained  by  scanning  all  data  channels  sequentially 
20  times  and  averaging  the  recorded  voltages  for  each  channel.  One 
scan  sequence  was  completed  in  approximately  200  ms.  The  average  volt- 
ages, percent  deviation  of  voltage,  and  reduced  data  were  displayed  on 
a cathode  ray  tube  terminal  for  operator  inspection.  The  display  was 
updated  every  5 seconds.  Selective  scans  were  saved  on  magnetic  tape 
and  in  a disc  file  for  prlnt-out  at  the  conclusion  of  the  run. 

Pitot  pressure  surveys  across  the  test  section  were  obtained  with 
an  analog  plotter  driven  by  two  digital-to-nnalcg  converters.  The  raw 
pressure  data  from  the  transducer  were  obtained,  digitized,  and  reduced 
to  pressure  ratio  (pitot  pressure/reservoii  pressure)  on-line  and  the 
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Figure  10.  Instrumentation  schematic 


results  were  displayed  on  the  plotter  in  real  time.  The  local  probe 
position  within  the  test  section  was  determined  from  the  voltages 
obtained  from  a potentiometer  connected  to  the  probe  drive  mechanisms 
and  were  used  to  drive  the  horizontal  axis  of  the  plotter. 

Wavelength  scans  were  obtained  by  displaying  the  digitized  photo- 
multiplier outputs  on  the  CRT  screen . The  spectrometer  was  scanned 
from  high  to  low  wavelengths  at  a constant  speed.  For  rotational  tem- 
perature measurements,  the  scan  speed  was  10  A/min  while  for  vibrational 
temperature  measurements  it  was  50  A/min.  The  corresponding  slit  widths 
were  50  and  150  m®>  respectively.  Reduced  rotational  temperatures  and 
integrated  vibrational  band  intensities  were  then  computed  on-line 
immediately  after  each  data  scan. 


IV.  RESULTS  AND  DISCUSSIONS 


A.  STATIC  CHAMBER  TESTS 

The  rotational  structure  of  the  bands  in  the  N^+  first  negative 
system  were  examined  in  detail  in  the  static  chamber  at  a static  tem- 
perature of  290  K and  with  pressures  ranging  from  0.2  to  25  torr.  Com- 
parative scans  of  the  rotational  structure  of  the  (0,0 ) band  of  the  N-_>+ 
first  negative  system  at  low  and  high  densities  are  shown  in  Fig.  11. 

The  only  apparent  effect  of  elevated  density  on  the  rotational  structure 
is  the  increased  intensity  of  the  rotational  lines  corresponding  to  K' 
values  of  16  19*  This  is  due  to  overlap  by  the  ( 3 ,6 ) band  of  the  N3 
second  positive  system.  Line-slope  plots  for  the  data  of  Fig.  11  are 
given  in  Fig.  12.  The  plot  of  tin  [I/Sk’  (g)v4  ] against  K'(K'  + 1)  should 
be  linear  with  a slope  inversely  proportional  to  the  rotational  temper- 
ature. The  low  density  data  display  such  linearity.  However,  the 
linearity  of  the  line-slope  plots  at  high  density  is  degraded  due,  in 
part,  to  the  overlap  by  the  (3,6)  second  positive  band. 

The  measured  rotational  temperatures  are  plotted  against  the  static 
pressure  in  Fig.  13.  The  indicated  rotational  temperature  Increases 
rapidly  from  the  known  value  (290  K)  as  the  pressure  is  increased, 
reaching  a nearly  constant  value  of  355  K at  pressures  over  3 torr.  'The 
scatter  of  the  data  in  Fig.  13  increases  rapidly  as  the  pressure  is 
raised  above  3 torr.  This  scatter  is  greater  than  that  attributed  to 
basic  uncertainties  in  determining  the  rotational  temperature  from  the 
line -slope  plots. 

Various  modifications  to  the  test  chamber  were  made  in  an  attempt 
to  reduce  the  rotational  temperatures  measured  at  high  densities.  In 
particular,  the  nitrogen  test  gas  was  injected  into  the  chamber  through 
a one-inch  diameter  copper  tube  oriented  with  its  axis  perpendicular  to 
the  electron  beam  and  with  the  tube  exit  within  approximately  0.5  in. 
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Figure  11.  Rotational  structure  of  the  N2+  first 
negative  (0,0)  band 
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Figure  12.  N2+  (0,0)  line-slope  plots  from  the  static  chamber 
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of  the  beam.  The  flow  velocity  through  the  tube  was  estimated  to  be 
approximately  30  ft/s.  This  was  done  in  an  attempt  to  remove  secondary 
electrons  and  low  energy  ions  created  by  the  beam  excitation  process 
frcm  the  region  of  observation  and  to  reduce  the  possibility  of  heating 
of  the  ambient  gas  by  the  electron  beam.  This  had  a negligible  effect 
on  the  measured  rotational  temperatures. 

Excitation  by  low  energy  particles  (secondary  electrons,  ions,  and 
metastable  molecules  and  atoms)  and  collision  quenching  are  the  most 
probable  causes  for  errors  in  rotational  temperature  at  high  density. 

Of  these  two,  excitation  by  low  energy  particles  is  the  most  likely. 

The  line-slope  plots  of  Fig.  12  display  reasonable  linearity  when  the 
overlap  of  the  (3,6)  band  of  the  second  positive  system  is  systematically 
removed  from  the  data.  In  addition,  visual  observations  of  the  gas 
within  the  test  chamber  at  high  densities  indicated  a large  degree  of 
excitation  occurred  in  regions  far  removed  from  the  electron  beam.  This 
visual  observation,  coupled  with  the  large  number  of  charged  particles 
collected  by  the  drift,  tube  of  the  electron  beam  generator  (which 
necessitated  control  of  the  potential  difference  between  the  drift  tube 
and  test  chamber),  support  the  conjecture  that  the  rotational  tempera- 
ture error  in  this  case  is  due  to  excitation  by  low  energy  particles. 

The  data  of  Ashkenas*  obtained  in  a static  chamber  at  a temperature 
of  300  K are  compared  with  the  present  data  in  Fig.  13.  The  correlation 
formula  developed  by  Ashkenas  over-estimates  the  rotational  temperature 
error  of  low  densities  and  is  in  qualitative  agreement  with  the  errors 
observed  here  at  high  densities.  However,  if  the  errors  are  due  to 
excitation  by  secondary  particles,  the  amount  of  excitation  will  depend 
strongly  on  the  geometry  of  the  test  chamber,  the  method  of  gas  injec- 
tion, the  electrical  condition  of  the  chamber  walls,  the  method  of 
collection  of  the  high  energy  electrons,  etc. 

It  is  notable  that  there  are  about  as  many  empirical  corrections 
to  convert  measured  rotational  temperatures  to  actual  gas  temperatures 
as  there  are  investigators  who  have  reported  such  results.  MHguire11 
obtained  different  errors  than  those  reported  by  Ashkenas'5  and  Kobben 
and  Talbot. *•'  Hickman  ’'5  modified  Muntz's  original  model  of  the  excitation- 
emission  process  in  an  attempt  to  describe  the  rotational  temperature 
errors,  particularly  for  temperatures  below  100  K.  However,  the  Hickman 
model  removes  only  some  of  the  error.  Maguire's  free  jet  data  indicate 
that  the  errors  are  due  to  selective  quenching  of  the  rotational  energy 
levels,  rather  than  excitation  by  low  energy  particles. 

The  cause  for  the  errors  in  rotational  temperature  measurement  at 
low  static  temperatures  is  still  unresolved  and  no  attempt  has  been  made 
here  to  contribute  to  the  understanding  of  this  problem.  However,  it  is 
clear  that  all  measurements  taken  in  a static  test  chamber  are  suspect 
due  to  excessive  excitation  of  the  N-/  first  negative  system  by  low 
energy  particles  which  remain  within  the  region  of  observation  for  un- 
determined lengths  of  time.  Since  the  energy  (velocity)  of  these 
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exciting  particles  is  .low,  the  excitation  process  occurs  over  a much 
greater  time  span  than  the  excitation  by  beam  electrons.  Hence,  there 
is  a high  probability  that  energy  will  be  exchanged  between  the  exciting 
particle  and  the  rotational  and  vibrational  energy  modes  of  the  molecule 
before  the  electronic  transition  occurs.  In  this  ease,  the  rotational 
and  vibrational  population  distributions  can  be  altered,  leading  to  in- 
dicated rotational  and  vibrational,  temperatures  which  bear  little  rela- 
tionship to  the  gas  transitional  temperature.  The  excitation  may  also 
result  in  a non -Boltzmann  population  distribution  so  that  a temperature 
for  the  affected  energy  mode  cannot  be  defined. 

The  ratio  of  intensities  of  the  (0,1)  and  (1,2)  bands  of  the  first 
negative  system  is  plotted  vs.  pressure  in  Fig.  14.  The  theoretical 
analysis  predicts  a band  intensity  ratio  equal  to  7.60.  However,  the 
ratios  measured  in  the  static  chamber  monatonically  increase  as  the 
pressure  is  raised. 

A dramatic  change  in  the  structure  of  the  (0,.l)  band  is  observed 
as  the  pressure  is  increased.  Comparing  the  low-density  profiles  of 
Fig.  2 with  the  same  band  intensity  distribution  at  20  torr  given  in 
Fig.  15  shows  a large  increase  in  the  apparent  intensity  of  the  R -branch. 
This  is  due  to  overlap  by  the  (.1,5)  band  of  the  N:?  second  positive 
system.  It  is  well  known  that  the  second  positive  system  is  excited 
primarily  by  low  energy  (3  eV)  electrons.11  Hence,  an  abnormally  high 
secondary  electron  concentration  trapped  within  the  static  test  chamber 
will  lead  to  large  intensities  in  the  second  positive  system. 

It  is  observed  (Fig.  15)  that  the  P-branch  of  the  (0,l)  band  is 
largely  unaffected  by  the  overlap  of  the  (.1,5)  band.  Hence,  the  ratios 
of  intensities  of  the  P-branch  portions  of  the  bands  were  investigated 
as  a means  of  determining  vibrational  temperatures  in  the  presence  of 
overlap  due  to  the  second  positive  bands.  Theoretical  predictions  of 
integrated  P-branch  intensity  ratios  were  made  with  results  similar  to 
those  given  In  Fig.  4. 

The  ratio  of  integrated  P-branch  intensities  for  the  (0,1)  and 
(.1,2)  bands  measured  in  the  static  test  chamber  are  shown  in  Fig.  16. 

Tile  P-branch  ratio  is  seen  so  be  nearly  constant  for  equivalent  pres  - 
sures  up  to  35  torr,  although  there  is  a large  amount  of  scatter  in  the 
data.  Ibis  scatter  is  due  primarily  to  difficulties  in  determining  the 
(.1,2)  P-branch  intensity.  The  I'-branch  intensity  was  assumed  to  include 
the  band  from  the  head  to  the  point  of  minimum  intensity  between  the  P 
and  R branches.  As  seen  in  Fig.  2,  the  division  between  the  two  branches 
is  well  defined  for  the  (0,1)  band  but  not  for  the  (1,2)  band.  However, 
since  the  integrated  P-branch  ratio  shows  no  trend  with  density,  it  may 
provide  a means  of  measuring  N-,  vibrational  temperature  in  the  presence 
of  band  overlap  providing  the  scatter  in  the  intensity  measurements  can 
be  reduced  to  an  acceptable  value. 
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Figure  15.  No+  first  negative  (0,1)  band  sequence  intensity 
measured  in  the  static  test  chamber 
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Data  obtained  by  Harbour'  are  compared  with  the  present  data  in 
Pigs,  lb  and  lu.  Harbour  performed  experiments  in  a static  chamber  with 
molecular  nitrogen  at  a temperature  of  296  K,  a beam  voltage  of  15  kV, 
and  beam  currents  near  0.1 5mA.  Both  total  band  intensities  and  band 

peak  intensities  were  measured  for  the  (0,1)  and  (1,2)  bands.  Harbour's 
data  display  the  same  basic  trend  as  reported  here  - an  increase  in 
intensity  ratio  with  increasing  gas  density.  Harbour's  ratio  of  band 
peaks  should  be  free  of  overlap  by  the  (i,5)  band  of  the  second  positive 
system  and  should  vary  with  gas  density  in  the  same  manner  as  does  the 
ratio  of  P-branch  intensities,  if  there  are  no  mechanisms  other  than 
band  overlap  which  contribute  to  the  variation  in  intensity  ratios  with 
density.  However,  Harbour's  ratio  of  band  peaks  increases  more  rapidly 
with  density  than  does  the  P-branch  intensity  ratio  of  the  present 
studies.  Since  the  rotational  temperature  measurements  show  an  apparent 
change  in  the  rotational  intensity  distribution  within  the  (0,0)  band  at 
high  densities,  it  is  likely  that  similar  changes  occur  in  all  bands  of 
the  first  negative  system.  Hence,  band  peak  intensities  and  integrated 
intensities  may  not  vary  in  the  same  fashion  as  the  gas  density  is  in- 
creased. 


It  is  concluded  that  the  static  test  chamber  experiments  are  of 
little  value  in  judging  the  overall  applicability  of  the  electron  beam 
method  at  high  density.  The  no-flow  conditions  within  the  chamber  lead 
to  excessively  high  concentrations  of  secondary  electrons,  low  energy 
ions,  and  molecules  in  metastable  energy  states.  These  particles  will 
excite  ground  state  nitrogen  molecules  with  a "collision-time"  which  is 
much  longer  than  that  associated  with  excitation  by  beam  electrons. 
During  the  relatively  long  excitation  time,  re-distribution  of  the 
energy  within  the  rotational  and  vibrational  energy  states  can  occur 
before  the  molecule  makes  a transition  to  the  N-»+B  E state.  Hence,  the 
basic  assumptions  of  the  excitation-emission  analysis  used  to  interpret 
the  observed  line  and  band  intensities  in  terras  of  rotational  and  vibra- 
tional temi>eratures  are  no  longer  applicable.  It  is  apparent  that 
reliable  calibration  experiments  can  be  conducted  only  in  flow  fields 
with  velocities  sufficient  to  sweep  the  low  energy  particles  out  of  the 
region  of  observation. 

In  the  present  studies,  the  apparently  large  contribution  of 
excitation  by  low  energy  particles  to  the  observed  radiation  which 
causes  significant  overlap  of  the  N?+  system  by  bands  of  the  N?  second 
positive  system  completely  overshadows  any  effects  of  selective  quench- 
ing of  the  particles  in  the  H2+B2I  state. 

B.  WIND  TUNNEL  STUDIES 

1.  Nominal  Run  Conditions 


For  nozzle  Nl,  wind  tunnel  data  were  obtained  at  nominal  axial 
locations  downstream  of  the  nozzle  exit  equal  to  0.10,  1.5,  and  2.0 
inches.  Since  much  lower  axial  gradients  in  flow  properties  existed  in 
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nozzle  N2,  data  were  collected  for  this  nozzle  only  at  n location  1.0 
inch  downs ti earn  of  the  nozzle  exit. 

The  1 . s and  2.0  inch  stations  of  N1  were  within  the  strong  expansion 
system  emanating  from  the  nozzle  exit  caused  by  the  large  pressure  un- 
balance between  the  nozzle  exit  and  test  cabin  static  pressures.  In  this 
region  of  the  flow  field,  large  axial  and  radial  gradients  in  all  flow 
properties  exist.  To  ensure  that  probe  and  electron  beam  data  were 
taken  at  the  same  position  within  the  flow  field,  the  electron  beam  was 
always  positioned  such  that  it  intersected  the  tip  of  the  probe  while 
the  probe  was  positioned  on  the  tunnel  centerline.  The  beam  position 
was  adjusted  with  a deflection  coil  mounted  between  the  exit  orifice  of 
the  drift  tube  and  the  test  region.  Once  the  proper  alignment  was 
achieved,  the  pitot  probe  was  withdrawn  and  spectral  data  Were  collected. 

To  evaluate  the  accuracy  of  the  electron  beam  measurements  taken  in 
the  arc-heated  facility,  a reliable  model  of  the  expansion  process  from 
the  reservoir  to  the  measuring  station  must  be  available.  Computations 
with  a computer  program  which  is  a modification  of  that  described  by 
Lordi26  showed  that  chemical  freezing  was  to  be  expected  early  in  the 
expansion  process.  Hence,  vibrational  relaxation  was  the  only  kinetic 
mechanism  of  importance. 

'theoretical  estimates  of  gas  properties  including  the  effects  of 
vibrational  relaxation  were  obtained  from  a modification  of  the  time 
dependent  method  of  Anderson  as  summarized  in  Ref,  27.  In  this  case,  the 
relaxation  of  the  vibrutional  energy  mode  is  assumed  to  be  described  by 
the  harmonic  oscillator  model.  That  is, 

dEv/dt  « (Eve-Ev)/Te  (26) 

where  Ev  is  the  local  vibrational  energy,  Eve  is  the  vibrational  energy 
evaluated  at  the  local  translational  temperature,  and  xe  is  the  effective 
vibrational  relaxation  time.  In  a multi-species  gas  it  is  common  to  evalu- 
ate xe  with  a "parallel  resistance"  mixture  rule.  For  example,  in  N^-0? 
mixtures,  xe  would  be  given  by 

lAe  = Xn?/t(N2-N2)  ♦ Xq2/t(N2-  02)  (27) 


In  mixtures  of  N2  and  02,  theoretical  and  experimental  results16  show 
that  both  gas  species  relax  with  the  same  effective  relaxation  time. 

Values  for  te  have  recently  been  determined  by  Sebacher16  for  N2 
relaxing  in  both  pure  N2  and  in  air.  The  N2-N2  and  N2-air  rates  are 
given  by 

PT(j2_N2  ::  3.4  x 10"lzexp(195T'l/3)  (28) 
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and 


pTN  -air  " l'2  x 10"10exp( 130T‘1/3 ) (29) 

where  x is  in  seconds,  p is  in  atmospheres,  and  T is  in  degrees  kelvin. 
These  rates  were  determined  in  Ref.  16  from  a set  of  measurements  ob- 
tained with  an  electron  beam  In  high  speed  N2  and  air  flows  at  gas  den- 
sities well  below  those  necessary  to  apply  the  electron  beam  method  with- 
out complications  due  tc  quenching  and  excitation  by  secondary  particles. 
Further  experiments  are  reported  in  Ref.  29  where  the  N?  rate  of  Ref.  16 
was  employed  with  the  computer  program  of  Ref.  27  to  establish  N2  vibra- 
tional relaxation  rates.  The  notable  feature  of  the  relaxation  rates 
given  in  Refs.  16  and  29  is  that  they  distinguish  between  relaxation  of 
N2  in  air  and  in  pure  N2.  'Ibis  is  contrasted  with  the  earlier  correla- 
tions of  Ref.  15  where  a single  vibrational  relaxation  rate  was  used  for 
both  gases.  The  rate  given  in  Ref.  15  to  be  applied  to  both  N2  and  air 
flows  is 

ip  = 3 x 10'1J  exp(l8l  T-1/3)  (30) 

The  relative  influences  of  the  rates  given  by  Eqs.  (28)  and  (30)  on 
the  predicted  N2  vibrational  and  gas  translational  temperatures  are  shown 
for  nozzle  N1  in  Fig.  17.  It  can  be  seen  that  the  single  rate  given  by 
Kq.  (30)  leads  to  less  vibrational  relaxation  and,  hcncc,  to  a higher 
vibrational  temperature  than  does  the  rate  of  Eq.  (28).  Since  different 
rates  of  N2  vibrational  relaxation  in  air  and  h'2  flows  are  anticipated, 
the  appropriate  rate  of  Eq.  (28)  or  (29)  was  used  in  all  calculation  of 
these  studies. 

The  computer  program  of  Ref.  27  was  employed  for  these  studies  be- 
cause it  includes  the  effects  of  vibrational  relaxation  on  the  entire 
flow  field;  that  is,  the  energy  changes  accompanying  vibrations],  relaxa- 
tion are  included  in  the  overall  energy  balance  for  the  system.  Ihe 
various  kinetic  rates  given  in  Ref.  28  were  employed  except  for  the 
N2-N2  rate.  The  rate  given  by  either  Eq.  (28)  or  (29)  was  used  in  place 
of  the  N2-N2  rate  and  the  vibrational  energy  of  O2  was  assumed  equal  to 
that  of  N2.  While  this  latter  assumption  gives  a poor  representation  of 
the  O2  vibrational  energy,  the  overall  effects  of  the  vibrational  re- 
laxation on  the  translational  energy  are  minor.  This  can  be  seen  in 
Fig.  17  where  the  effects  of  varying  the  N ? rate  lead  to  very  little 
change  In  the  translational  temperature. 

The  conditions  of  Fig.  17  correspond  to  the  highest  reservoir 
temperature  encountered  before  significant  oxygen  dissociation  occurs  in 
the  flow  field.  As  the  reservoir  temperature  decreases,  the  amount  of 
vibrational  relaxation  present  in  the  flow  decreases,  as  does  the  sensi- 
tivity to  the  particular  rate  equation  used.  For  many  of  the  wind 
tunnel  run  conditions  of  these  studies,  little  differences  exist  between 
the  test-section  vibrational,  and  translational  temperatures  predicted 
with  any  of  the  rate  equations. 


40 


NOZZLE  CONTOUR 


At  reservoir  temperatures  much  higher  than  that  of  Fig.  17,  signifi- 
cant oxygen  dissociation  occurs  so  that  the  assumption  of  equal  nitrogen 
ar.d  oxygen  vibrational  energies  has  even  less  effect  on  the  overall 
energy  balance.  The  main  effect  of  high  reservoir  temperatures  is  to 
decrease  the  oxygen  concentration  so  that  the  second  term  in  Eq.  (27) 
becomes  negligible.  In  this  case,  the  rate  given  by  Eq.  (28),  rather 
than  that  of  Eq.  (20),  should  be  used. 

To  implement  the  computer  program  of  Ref.  27  for  the  present  studies, 
the  CO2  and  H20  concentrations  inherent  in  the  program  must  be  set  to 
small  values.  Mo  difficulties  were  encountered  in  setting  the  H20  mole 
fraction  to  zero.  However,  numerical  difficulties  were  encountered  when 
the  CO2  concentration  was  set  Lo  zero.  A brief  set  of  numerical  experi- 
ments was  performed  to  establish  a minimum  mole  fraction  of  CO2  which 
would  allow  execution  of  the  program.  This  value  was  determined  to  be 
approximately  10-15.  Increases  in  the  minimum  value  by  at  least  seven 
orders  of  magnitude  had  no  effect  on  the  N2  vibrational  and  translational 
temperatures  to  four  significant  digits. 

When  the  low  area  ratio  wind  tunnel  nozzle  (Nl)  was  used,  the  high 
pumping  capacity  of  the  vacuum  system  led  to  a strong  expansion  system 
at  the  nozzle  exit  for  all  run  conditions.  Since  the  expansion  system 
causes  very  rapid  decreases  in  flow  temperature  and  density,  it  is  un- 
likely that  the  vibrational  energy  mode  will  relax  sufficiently  fast 
enough  to  respond  to  these  rapid  changes  in  flow  properties . Hence,  the 
vibrational  relaxation  was  assumed  to  proceed  only  downstream  to  the 
point  where  the  leading  expansion  wave  from  the  nozzle  exit  crossed  the 
nozzle  centerline.  All  vibrational  temperatures  downstream  of  this 
freezing  point  were  assumed  to  be  equal  to  their  values  at  the  freezing 
point.  As  can  be  seen  in  Fig.  17,  the  effect  of  this  freezing  is  rela- 
tively minor  and  is  important  only  for  the  data  obtained  at  the  farthest 
downstream  axial  position.  For  the  high  area  ratio  nozzle  (N2)  vibra- 
tional relaxation  was  assumed  to  proceed  from  the  reservoir  to  the 
measuring  station. 

The  gas  conditions  at  the  test  stations  were  determined  by  measur- 
ing the  local  pitot  pressures  (pT2).  The  theoretical  model  of  Ref.  27 
was  modified  to  include  a pitot  pressure  calculation.  This  was  accomp- 
lished by  assuming  that  the  flow  through  the  pitot  pretoe  shock  wave  was 
vibrationally  frozen.  A frozen  speed  of  sound  was  calculated  from 

af  = (YR/MT)1/2  (31) 

where  y is  the  frozen  ratio  of  specific  heats,  M is  the  molecular  weight 
of  the  gas  mixture,  and  T is  the  local  static  temperature.  The  Mach  number 
obtained  fran  af  and  the  local  velocity  was  then  used  in  the  isentropic 
formula  for  pTz/p\,  where  pj  is  the  static  pressure  upstream  of  the  probe 
shock  wave.  The  ratio  of  pitot  and  reservoir  pressures  was  then  obtained 
from  pT2/p0  **  pT2/pj  x Pj/Po?  where  p0  is  the  reservoir  pressure.  The 
pitot  pressure  ratios  measured  on  the  tunnel  centerline  were  matched  with 
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those  of  this  theoretical  model  to  give  all  flow  conditions  at  the 
measuring  stations. 

Typical  pitot  pressure  surveys  taken  at  the  0.1  and  2.0  inch 
stations  of  nozzle  N1  are  given  in  Fig.  18a  while  a typical  survey  for 
nozzle  N2  is  given  in  Fig.  18b.  Representative  reservoir  conditions  and 
predicted  test  section  gas  properties  are  summarized  in  Table  II. 

2.  Overall  Spectral  Features 

A wavelength  scan  of  the  entire  spectral  region  of  interest  is 
shown  in  Fig.  19.  These  data  were  obtained  at  the  0.10  inch  axial 
station  of  N1  at  a reservoir  temperature  of  11)00  K,  a reservoir  pressure 
of  0.8  atm,  and  an  equivalent  pressure  near  20  torr.  Hence,  they  corres- 
pond to  the  maximum  gas  density  investigated  in  these  studies.  The  scan 
is  uncorrected  for  the  spectral  response  of  the  optical  system.  Vibra- 
tional bands  due  to  the  N2+  first  negative,  N2  second  positive,  and  NO  y 
systems  dominate  the  spectrum.  Bands  from  the  0?+  second  negative  system 
were  also  observed  but  with  intensities  much  less  than  those  of  the  N2 
and  NO  bands . 

The  radiation  intensities  observed  at  high  densities  are  compared 
with  those  at  the  lower  densities  obtained  with  nozzle  N2  in  Fig.  20. 

The  most  notable  difference  between  these  two  scans  is  the  obvious  appear- 
ance of  relatively  intense  radiation  in  the  02+  second  negative  bands 
(particularly  the  (3»3)  band)  in  the  lower  density  data.  This  implies  a 
significant  difference  in  the  relative  quenching  of  the  NO y and  02+  second 
negative  bands  at  high  density  which  leads  to  severe  overlap  of  the  02+ 
bands  by  those  of  the  NO  y system.  Sensitivity  experiments  were  conducted 
to  determine  the  operating  conditions  which  allow  02  vibrational  tempera- 
ture measurements  to  be  made  in  the  presence  of  NO  and  are  discussed  below. 

3.  N2  Results 

The  rotational  structures  of  the  (0,0)  and  (0,1)  bands  of  the  N2+ 
first  negative  system  were  examined  In  detail  to  determine  if  excitation 
by  secondary  particles  and  preferential  collision  quenching  disturb  the 
band  structure  in  a flow  with  high  gas  velocity.  Scans  of  the  R-branches 
of  the  (0,0)  and  (0,1)  bands  obtained  in  nozzle  N1  at  an  equivalent 
pressure  of  approximately  20  torr  and  a rotational  temperature  of  786  K 
are  show/!  Figs.  21  and  22.  The  changes  in  the  intensities  of  the 
rotations,  lines  in  the  (0,0)  band  with  rotational  quantum  numbers  greater 
than  15  due  to  overlap  by  the  (3,6)  band  of  the  N2  second  positive  system 
are  barely  discernible  in  Fig.  21.  In  addition,  there  is  no  evidence  of 
significant  overlap  of  the  (0,1)  band  by  the  (1,5)  band  of  the  second 
positive  system. 

The  line-slope  plot  for  the  data  of  Fig.  21  is  shown  in  Fig.  23.  The 
departure  from  linearity  is  due  primarily  to  interference  of  the  P-branch 
lines  with  the  R-branch  lines  caused  by  the  high  rotational  temperature. 
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PITOT  PRESSURE  (PSIA) 


i 


Figure  18a.  Typical  pitot  pressure  surveys  in  nozzle  Nl. 

X = distance  downstream  of  nozzle  exit  in  inches. 
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Figure  18b.  Typical  pitot  pressure  survey  in  nozzle  N2; 
1.0  inch  downstream  of  nozzle  exit. 


Table  II.  Nominal  Run  Conditions 


Inches  Downstream 
Nozzle  Exit 

To 

(K) 

Po 

(atm) 

P 

(torr) 

PEQIV 

(t.orr) 

Mach 

No. 

Nozzle  Nl: 

0.1 

2100 

0.66 

38.04 

1006 

U.3 

2-3 

1700 

O.95 

59-3 

809 

22.0 

2.3 

1.5 

1500 

1.02 

14.0 

476 

8.8 

3-3 

2.0 

2600 

0.82 

4.7 

648 

2.2 

4.0 

2100 

0.70 

1.63 

401 

1.2 

4.6 

1700 

O.89 

i.ge 

318 

1.8 

4.7 

Nozzle  N2: 

1.0 

550C 

1.20 

0.234 

579 

0.121 

6.7 

3700 

0.74 

0.144 

381 

0.113 

6.8 

2900 

0.60 

0.113 

284 

0.119 

6.8 
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INTENSITY  INTENSITY 


(3,3)  02+(2-) 

I 


WAVELENGTH  (A) 


a.  Nozzle  Nl;  Pgg  « 20  torr 


WAVELENGTH  (&) 


b.  Nozzle  N2;  PgQ  * 0.18  torr 


Figure  20.  Spectral  scans  of  electron  beam  induced  emission 
in  the  arc  heated  wind  tunnel. 
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INTENSITY 


Figure  21.  No+  (0,0)  band  rotational  structure  from  the 
wind  tunnel;  peq  = 20  torr;  nozzle  Nl. 
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INTENSITY 


WAVELENGTH  (A) 


Figure  22.  N,-;  (0,1)  band  rotational  structure  frcta  the 

wind  tunnel;  pgQ  = 20  torr;  nozzle  Nl. 


Figure  2J.  Line  slope  plot  for  data  of  Figure  21 . 
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However,  lines  in  the  tail  of  the  band  which  are  relatively  free  of 
interference  from  the  P-branch  lines  show  reasonable  linearity. 

Analyses  of  these  spectral  details  show  that  the  excitation  by 
secondary  particles  which  dominate  the  measurements  in  the  static 
chamber  at  these  densities  are  not  observable  in  the  wind  tunnel  data. 

In  addition,  the  significant  overlap  by  the  bands  of  the  second  positive 
system,  which  are  excited  primarily  by  secondary  electrons,  do  not 
appear  to  contaminate  the  first  negative  bands  in  any  significant  manner. 

A summary  of  the  N2  data  obtained  from  the  wind  tunnel  experiments 
is  given  in  Table  III.  The  measured  vibrational  temperatures  are  compared 
with  the  theoretical  predictions  in  Fig.  24.  The  maximum  scatter  in  the 
data  of  Fig.  24  amounts  to  _+  61  and  is  due  to  accumulation  of  errors  in 
determining  the  reservoir  temperature  from  the  measured  mass  flow  rates 
and  reservoir  pressure,  uncertainties  in  determining  the  rotational 
temperature  from  the  line  slope  plot  (for  use  with  Fig.  I to  determine 
vibrational  temperature),  and  inaccuracies  in  obtaining  integrated  band 
intensities . 

A 51  error  in  mass  flow  rate  will  result  in  a 10^  error  in  reservoir 
temperature.  For  these  studies,  errors  in  the  measurement  of  both  the 
primary  and  secondary  mass  flow  rates  contribute  to  the  total  error  in 
reservoir  temperature.  Errors  in  determining  the  rotational  temperature 
from  the  line  slope  plot  are  due  primarily  to  uncertainties  in  fitting 
a straight  line  through  the  data.  At  high  rotational  temperatures 
(800  K),  the  interference  of  P-branch  lines  with  the  R-branch  lines  makes 
this  determination  more  difficult  than  it  is  at  lower  temperatures. 

Under  these  conditions,  the  accuracy  in  measured  rotational  teirperature 
Is  no  better  than  approximately  ♦_  6JJ.  Hence,  the  scatter  of  the  data 
evident  in  Fig.  24  is  within  the  combined  accuracy  with  which  the 
measured  and  predicted  vibrational  temperatures  can  be  determined  with 
the  present  system. 

The  coiqparisons  of  the  theoretical  and  experimental  N2  vibrational 
temperatures  shown  in  Fig.  24  demonstrate  the  effect  of  gas  composition 
on  the  rate  of  N2  vibrational  relaxation.  In  air  flows,  the  te  given  by 
Eq.  (27)  describes  the  rate  of  relaxation  of  both  the  N2  and  O2.  Its 
value  is  less  than  that  for  pure  M2  but  greater  than  that  for  pure  02 . 

As  the  O2  concentration  decreases  due  to  dissociation,  xe  increases, 
leading  to  less  vibrational  relaxation  in  the  flow  field. 

The  02  mole  fractions  at  the  various  test  conditions  of  these 
studies  are  shown  as  functions  of  the  reservoir  temperatures  in  Fig.  24. 
The  data  obtained  with  no/szle  N1  correspond  to  O2  mole  fractions  near 
0.20,  typical  of  low  temperature  air  and  the  N^-air  rate  of  Eq.  (29)  is 
expected  to  apply.  However,  the  test  conditions  of  nozzle  N2  resulted 
in  significant  oxygen  dissociation  so  that  the  N2-N2  rate  of  2q.  (28) 
should  apply,  parti cularly  at  the  higher  reservoir  temperatures.  The 
measured  N2  vibrational  temperatures  shown  in  Fig.  24  agree  with  these 
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Table  III.  Run  Summary  - Ng  Data 


Rvui 

To 

00 

Po 

(atm) 

P 

(torr ) 

PEQ 

(torr) 

Tr(Nj) 

(K) 

tv(n-,) 

(K) 

tv(thbor 

(K) 

Nozzle 

Nl: 

52 

2561 

0.799 

5.00 

2.18 

600 

2182 

2135 

60 

2106 

O.658 

39.5 

11.3 

902 

1897 

1914 

61 

1814 

1.047 

61.6 

20.7 

839 

1634 

1647 

64 

2104 

0.724 

43.6 

12.5 

- 

1740 

I898 

68 

2136 

0.704 

1.68 

1.19 

453 

2000 

1924 

1766 

o.89e 

2.03 

1.78 

371 

1701 

1637 

70 

1456 

1.020 

14.48 

9.18 

626 

1358 

1377 

85 

2086 

0.688 

4.30 

2.62 

49« 

1690 

1390 

86 

1497 

1.024 

6.40 

4.81 

398 

1430 

1416 

2121 

0.671 

4.20 

2.52 

498 

1780 

1909 

87 

1763 

0.652 

38,4 

13.4 

915 

1620 

1663 

1392 

1.024 

59.0 

26.1 

754 

1189 

1331 

88 

1973 

0.672 

2.23 

1.56 

452 

1880 

1817 

1616 

1.048 

3.38 

2.42 

419 

1590 

1510 

Nozzle 

N2: 

3 

4258 

1.032 

0.267 

0.179 

449 

3300 

3500 

5 

2866 

0.594 

0.111 

0.119 

279 

2750 

2650 

2974 

0.615 

0.116 

0.119 

292 

3000 

2730 

6 

3764 

0.714 

0.139 

0.107 

390 

3250 

3252 

3685 

0.700 

0.136 

0.107 

380 

3430 

3225 
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Table  III.  Run  Summary  - N?  Data  - (Continued) 


Run 

To 

(K) 

Po 

( atm) 

P 

( torr ) 

PEQ 

( torr ) 

tr(n2) 

(K) 

Tv( N2 ) 
(K) 

Tv(  THEORY ) 
(K) 

Nozzle  N2:  (Continued) 

7 

3680 

0.744 

0.144 

0.114 

379 

3260 

3220 

3725 

0.756 

0.147 

0.115 

384 

3330 

3250 

8 

3686 

0.988 

0.192 

0.152 

380 

3230 

3230 

3825 

1.013 

0.198 

0.149 

397 

3300 

3280 

9 

3592 

0.908 

0.176 

0.143 

368 

3120 

3150 

3636 

0.913 

0.177 

0.142 

373 

26 0C 

3180 

10 

3427 

0.745 

0.143 

0.124 

346 

3270 

3005 

3514 

0.768 

0.148 

0.124 

358 

3240 

3100 

12 

3336 

0.760 

0.145 

0.130 

336 

3070 

3000 

13 

3541 

1.036 

0.200 

0.166 

361 

3110 

3110 

3584 

1.061 

0.205 

0.168 

367 

3220 

3130 
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Figure  24.  Measured  and  predicted  Nj  vibrational  temperature; 

O nozzle  Nl;  Q nozzle  N2;  —theoretical  Tv; 

— - theoretical  t 02- 


05 


Irent:.  The  low  tempera Lure  data  agree  with  theoretical  prediction;: 
employing  tlie  Nj-air  rate  while  the  high  temperature  data  are  in 
better  agreement  w i t.h  theoretical  results  obtained  wl  t.h  the  N;--M7  rate. 

The  analysis  of  the  N-.  vibrational  relaxation  process  used  here 
allows  for  only  one  effective  collision  partner  in  de-activating  the  N? 
vibrational  energy  mode.  A more  general  analysis  should  be  formulated 
employing  a rate  in  the  form  of  Eq . (2?)  which  allows  for  varying  0;> 
concentrations  throughout  the  flow  field.  Howevei , the  aim  of  the  pre- 
sent studies  is  to  determine  the  accuracy  with  which  vibrational  tempera- 
tures can  be  measured  over  a wide  range  of  flow  densities.  Hence,  it 
was  considered  sufficient  to  obtain  theoretical  estimates  of  f^ow  pro- 
perties with  the  N,> -air  rate  at  the  low  reservoir  temperatures  of  nozzle 
N1  while  the  No-No  rate  was  used  for  the  high  reservoir  temperatures  of 
nozzle  N2. 

To  investigate  the  possible  effects  of  gas  density  on  the  measured 
N;>  vibrational  temperatures,  the  ratios  of  the  measured  and  predicted 
vibrational  temperatures  listed  in  Table  III  are  plotted  in  Fig.  25 
against,  equivalent,  pressure.  No  specific  trend  in  the  ratio  of  the 
measured  and  predicted  temperatures  with  gas  density  is  evident.  At 
equivalent  pressures  greater  than  approximately  1 torr,  the  radiation  in 
the  tl,*  first  negative  system  is  highly  quenched  by  colltsional  mechan- 
isms. However,  the  data  of  Fig.  25  demonstrate  that  the  rotational  and 
vibrational  energy  levels  in  the  N;>*(B?E)  electronic  energy  state  are 
not  preferentially  quenched  nor  are  i.m,-  data  influenced  by  excitation 
by  secondary  particles  in  any  measurable  amount  for  equivalent  pressures 
at  least  up  to  20  torr. 


A.  Nitrio  Oxide  Results 


Ttie  ratios  of  intensities  of  the  NO  y (0,4)  and  (1,5)  bands 
obtained  in  the  wind  tunnel  experiments  ure  summarized  in  Table  IV.  The 
band  Intensity  ratios  are  scaled  by  the  value  at  300  K as  discussed  pre- 
viously to  account  for  variations  in  the  spectral  s<..isl tivity  of  the 
optical  system.  The  corresponding  Nj  theoretical  vibrational  tempera- 
tures are  included  to  give  an  indication  of  the  NO  vibrational  tempera- 
ture. The  theoretical  N?  vibrational  temperature*;  are  used  because  of 
the  uncertainty  In  the  vibrational  relaxation  rate  for  NO  in  expansion 
flows.  While  experimental  excitation  rates  for  NO  are  available  from 
shock  tube  studies,  the  applicability  of  these  rates  in  expansion  flows 
is  uncertain.  It  is  well-known13-16  that  vibrational  excitation  rates 
for  N2  greatly  underestimate  the  vibrational  relaxation  in  expansion 
flows.  It  is  likely  that  similar  effects  occur  in  expanding  NO. 

The  (0,4)  and  (1,5)  band  intensity  ratios  obtained  at  all  run 
conditions  are  plotted  in  Fig.  26  against  the  N.,  tlieoretical  vibrational 
temperature.  Included  In  Fig.  26  are  the  theoretical  band  intensity 
ratios.  The  theoretical  band  intensity  ratio  for  direct  excitation  is 
"Iso  shown  scaled  such  that  Its  mom  temperature  value  is  equal  to  that 
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EQUIVALENT  PRESSURE.  3RR) 

Figure  25.  Ratios  cf  measured  and  predicted  N2  vibrational 
temperature;  O nozzle  Nl;  grange  of  data  for 
nozzle  N2  (except  second  condition  of  Run  9) 


Table  IV.  Run  Nummary  - Nitric  Oxide  Data 


Run 

To 

(K) 

Po 

(atm) 

P 

(torr ) 

PEQ, 

(torr) 

Tr(No) 

(K) 

tv(n2>  theory) 
(K) 

Wlis 

Nozzle  Nl: 

73 

1488 

0.942 

19.6 

10.9 

537 

l4l6 

1.918 

77 

2083 

0.756 

16.4 

6.4 

776 

1880 

1.847 

79 

1561 

O.983 

74.1 

27.3 

815 

1472 

1.989 

8l 

1654 

0.997 

2.4l 

2.22 

325 

1547 

2.088 

85 

2063 

O.687 

4.29 

2.64 

487 

1870 

2.068 

86 

1572 

l.OUl 

6.51 

4.66 

418 

1487 

2.119 

2121 

0.670 

4.20 

2.52 

500 

1909 

1.964 

87 

1.807 

0.671 

39-5 

13.5 

939 

1704 

1.645 

1392 

1.024 

59-0 

26.1 

754 

1331 

1.835 

88 

1973 

0.672 

2.23 

1.56 

452 

1817 

1.700 

I6l6 

1.048 

3.38 

2.42 

419 

1510 

2.035 

Nozzle  N2: 

3 

4053 

1.012 

0.199 

0.140 

425 

3400 

1.869 

3155 

1.019 

0.138 

0.138 

437 

3450 

1.681 

5 

2949 

0.609 

0.114 

0.119 

289 

2730 

2.052 

2986 

0.614 

0.115 

0,118 

294 

2750 

2.064 

6 

3679 

0.704 

0.137 

0.108 

380 

3240 

2.027 

3770 

0-707 

0.137 

0.106 

388 

3280 

2.108 

7 

3680 

0-753 

0.146 

0.116 

379 

3220 

1.968 

3725 

0-755 

0.147 

0.115 

384 

3240 

1.752 
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Table  IV.  Run  Summary  - Nitric  Oxide  Data  (Continued) 


Run 

To 

(K) 

Po 

(aim) 

P 

( torr ) 

PEQ 
( torr ) 

Tr(N2) 

(K) 

TV(N2, THEORY) 

(K)  WllS 

Nozzle  N2:  (Continued) 

8 

3732 

0.982 

0.191 

0.149 

385 

3230 

1.766 

3778 

0.992 

0.193 

0.148 

391 

3260 

1.732 

9 

3592 

0.910 

0.176 

0.143 

368 

3160 

1.953 

3592 

0.911 

0.176 

0.143 

368 

3160 

1.917 

10 

3427 

0.746 

0.143 

0.124 

347 

3050 

2.010 

3427 

0.746 

0.143 

0.124 

347 

3050 

2.030 

3514 

0.769 

0.148 

0.124 

358 

3100 

1.939 

12 

3336 

0.760 

0.145 

0.129 

339 

3000 

2.173 

3302 

0.756 

0.144 

0.130 

332 

2990 

2.180 

13 

3541 

1.052 

0.203 

0.169 

360 

3110 

2.046 

3585 

1.057 

0.204 

0.167 

366 

3140 

1.933 
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(0,41/0,5)  INTENSITY  RATIO 


Figure  26.  Comparison  of  measured  and  predicted  band  profiles 
for  NO  y (0,4)  and  (1,5)  bands; 

direct  excitation;  vibrational  relaxation 

model;  O nozzle  Nl;  □ nozzle  N2. 


from  the  vibrational  relaxation  model.  This  scaling  has  been  discussed 
previously  with  reference  to  Fig.  7 and  is  equivalent  to  non- dime  ns  ion- 
aliziing  the  measured  intensities  by  their  room- temperature  values. 

The  direct  excitation  model  (unsealed)  greatly  under-estimates  the 
band  intensity  ratios  at  all  run  conditions.  Hence,  the  effects  ob- 
served in  the  static,  chamber  are  repeated  even  when  high  gas  velocities 
are  present  within  the  beam  excitation  region. 

The  band  intensity  ratios  scaled  by  their  room  temperature  values 
obtained  with  nozzle  N1  under  high  density  conditions  are  under-estimated 
by  the  vibrational  relaxation  theory.  However,  the  agreement  between  the 
theoretical  and  experimental  ratios  at  the  lower  density  conditions  of 
nozzle  N2  is  much  better.  In  this  case,  the  direct  excitation  model 
gives  as  good  a representation  of  the  scaled  ratios  as  does  the  vibra- 
tional relaxation  model. 


The  percent  difference  between  the  theoretical  an  experimental 
band  intensity  ratios  are  plotted  in  Fig.  27  against  the  rotational 
temperature.  It  has  been  assumed  that  the  N2  and  NO  rotational 
temperatures  are  equilibrated  and  are  equal  to  the  tranlational 
temperature.  Except  for  a few  isolated  points,  a definite  trend  in 
the  difference  is  evident  in  Fig.  27.  In  addition,  there  is  no  appar- 
ent effect  of  the  gas  density  level  on  the  differences  between  the 
theoretical  and  experimental  ratios  of  band  intensities  shown  in  Fig, 

27 , The  correlation  of  the  difference  with  rotational  temperature  in- 
dicates that  the  theoretical  model  does  not  give  the  correct  intensity 
distribution  for  the  rotational  lines  within  the  bands.  This  can  be 
observed  more  clearly  in  the  comparative  band  profiles  of  Fig.  23. 

The  most  notable  difference  between  the  experimental  profile  and 
that  obtained  with  the  relaxation  model  shown  in  Fig.  28  is  that  the 
experimental  results  display  a larger  amount  of  energy  at  high  rotational 
quantum  numbers  in  both  bands.  This  energy  is  sufficient  to  introduce 
measurable  overlap  of  the  (0,4)  band  with  the  (1,5)  band.  The  inappli- 
cability of  the  direct  excitation-emission  model  is  also  evident  in 
Fig.  28. 

The  difference  between  the  theoretical  and  experimental  band  pro- 
files are  not  large  considering  the  assumptions  employed  in  the  theory. 
The  transition  probability  (C]j/Ca)kio  is  determined  from  room-tempera- 
ture, static  calibrations  and  the  results  are  applied  to  flows  with 
varying  rotational  and  vibrational  temperatures.  As  seen  in  Fig.  28, 
the  model  appears  to  give  a good  representation  of  the  relative  band 
intensities.  However,  the  failure  of  the  model  to  adequately  describe 
the  distribution  of  intensities  within  the  bands  is  symptomatic  of  the 
omission  of  an  important  excitation  or  de-excitation  mechanism. 

Since  the  y bands  result  from  resonant  transitions,  self-absorption 
of  the  emitted  radiation  may  provide  an  important  means  of  excitation. 
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Figure  27.  Percent  difference  in  theoretical  and  experimental  intensity  ratios  for 
the  NO  y (0,4)/(i,5)  bands;  O nozzle  Nl;  Q nozzle  N2. 
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Figure  28.  Comparisons  of  measured  and  predicted  band 
profiles  for  NO  y (0,4)  and  (1,5)  bands. 
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The  studies  cf  Ref.  5 show  for  bands  other  than  those  involving  the 
ground  vibrational  energy  level  in  the  2n  state  that,  self-absorption  will 
not  alter  the  emitted  radiation  intensities.  However,  the  cumulative 
effects  of  self-absorption  in  providing  a path  for  excitation  to  the  2I+ 
may  be  significant.  This  effect  can  be  modeled  in  a crude  way  by  expres- 
sing the  rate  of  excitation  to  a particular  vibrational  energy  level  in  the 
2tt  state  as 

<t>(v')  = ve  )Xnq  + hev3An*R(n*,X{g0)  (31) 

where  n*  is  the  number  density  of  excited  species  in  the  v'  level  of  the 
2>;  state,  R(n*,X^0)  is  the  rate  of  reabsorption  by  ground  state  NO  of 
photons  emitted  in  fluorescence,  and  all  other  terms  have  meanings  as 
previously  defined.  The  radiation  intensity  resulting  from  this  excita- 
tion can  be  written  as 


y-  = 1 - R(n*,XNQ)/n*A 

where  I0  is  the  intensity  in  the  absence  of  absorption. 


(32) 


Melton  and  Klemperer30  show  that 


R(n*,XN.0)/n*A  = 


OiPwqV?)1,2 

0.33  + (kiP^t/2)1-2 


(33) 


Eq.  (33)  is  formulated  for  population  of  a single  rotational  energy  level 
and  the  index  i runs  over  all  branches  in  the  band  originating  from  the 
level.  (Aj/A)  is  the  fraction  of  photons  emitted  in  the  i-th  branch  of 
the  ( v ' , 0 ) band . 


The  analysis  leading  to  Eq.  (33)  applies  only  when  the  gas  can  be 
considered  as  an  infinite  slab  of  thickness  t with  atomic  line  radiation 
(rather  than  band  radiation)  normal  to  the  slab.  The  quantity  kj  is  the 
absorption  coefficient  of  the  i-th  branch  and  Pmq  is  the  NO  partial  pres- 
sure. Ref.  30  gives  estimates  for  R(n*,Xf$Q)  which  show  that  in  pure  NO 
at.  a few  tenths  of  a torr  the  intensity  can  be  increased  by  as  much  as 
60%  by  self-absorption. 


Inclusion  of  self -absorption  in  the  analytical  model  of  tho  excita- 
tion-emission process  requires  terms  similar  to  that  of  Eq.  (33)  for  each 
rotational  line  in  the  NO  spectrum.  Further,  a reulistlc  representation 
of  the  intensity  distribution  within  the  electron  beam  would  be  required 
to  adequately  describe  the  absorption  process.  If  such  a theoretical 
model  could  be  developed,  it  might  not  only  describe  the  intensity  dis- 
tributions within  the  bands  but  might  also  eliminate  the  need  to  include 
upper-state  vibrational  relaxation  to  account  for  the  observed  relative 
band  intensities. 
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5.  Molecular  Oxygen  Results 

Comparative  scan  of  the  spectral  region  near  the  0?+(3>3)  and 
(4,3)  bands  obtained  from  the  wind  tunnel  and  static  chamber  tests  are 
shown  in  Fig.  29.  The  wind  tunnel  data  from  nozzles  N1  and  N2  were 
obtained  at  equivalent  pressures  of  4.20  torr  and  0.14  torr,  respectively. 
The  static  chamber  data  were  obtained  at  a pressure  of  0.40  torr  in  pure 
02* 


At  the  high  density  conditions  of  nozzle  Nl,  the  O^^C  3,3)  band  is 
completely  overlapped  by  the  NO  y (0,5)  and  N2  second  positive  (3,0) 
bands.  However,  at  the  low  density  conditions  of  nozzle  N2,  the  overlap 
is  much  less  seveie.  In  fact,  in  the  region  near  the  02+(3>3)  band  the 
low  density  wind  tunnel  data  resemble  very  closely  that  obtained  in  the 
static  chamber  with  pure  02. 

The  variation  of  the  overlap  of  the  02 + bands  with  changing  density 
indicates  that  the  02+  radiation  is  quenched  more  effectively  than  is  the 
NO  y radiation  as  the  density  increases.  To  investigate  this  effect  in 
more  detail,  a series  of  sensitivity  experiments  was  conducted  in  the 
static  test  chamber.  It  was  determined  that  the  minimum  beam  current- 
equivalent  0?  pressure  product  for  measurable  intensities  in  the  02  (4,3) 
band  is  approximately  0.1  torr-mA,  ignoring  any  effects  of  collision 
quenching. 

The  (4,3)  band  is  observable  only  in  the  static  chamber  data.  In 
the  high  density  data,  it  is  completely  overlapped  by  the  tall  of  the 
NO  y (1,6)  band.  However,  for  the  high  density  wind  tunnel  conditions  of 
Fig.  29,  the  current-02  pressure  product  was  0.26  torr-mA.  Since  tills 
vulue  is  well  ab(  ve  the  quenching  free  sensitivity  level  of  0.1  torr-mA, 
the  ubsence  of  detectable  02+  radiation  in  the  (4,3)  band  indicates  that 
the  02+  system  is  highly  quenched  at  the  higher  density  levels.  In  fact, 
the  NO  y band;  arc  so  much  more  intense  than  the  02  bands  that  it  is 
difficult  to  separate  the  02+  bands  from  the  overlapping  NO  y radiation. 

"Vie  02+( 4,3)  band  also  could  not  be  observed  in  the  low  density 
wind  funnel  data.  However,  the  O2  concentration  in  the  low  density  data 
was  approximately  1.6  x 10*4  cm-3  (Po2  = .007  torr)  giving  a curront-02 
pressure  product  of  approximately  .07  torr-mA;  hence,  the  absence  of 
observable  radiation  in  the  0?+(4,3)  band  in  these  wind  tunnel  experiments 
agrees  with  the  results  of  the  sensitivity  analyses. 

The  radiation  in  the  02 ! second  negative  bands  at  the  low  density 
wind  tunnel  conditions  was  clearly  distinguishable  from  the  other  band 
systems  present  in  the  spectrum.  In  addition,  the  radiation  intensities 
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Figure  29.  Spectral  scans  for  wind  tunnel  and  static  chamber  conditions . 
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agreed  will)  those  reported  earlier7’8.  At  the  high  density  conditions 
of  no/.?, le  N.L  ( Pj'Q  > 2.0  t.orr),  the  severe  overlap  of  the  07+  band’s  by 
the  NO  y bands  precludes  use  of  the  0->+  system  for  measuring  02  tempera- 
tures and  densities  with  the  spectral  resolution  of  these  studies  (0.2a 
max)  and  electron  hean  currents  at  Least  up  to  10  niA. 


C . UKAM  SHAPE 

An  important  concern  in  utilizing  an  electron  beam  is  the  spatial 
extent  of  the  luminescence.  For  density  measurements  it  is  imperative 
that  the  entire  beam  be  viewed  at  a given  location  along  the  beam  so 
that  calibrations  of  intensity  vs.  density  at  a Known  beam  current  can 
be  applied.  As  the  beam  spreads  with  increasing  density  and  beam 
length,  the  spatial  resolution  of  density  measurements  will  degrade. 

For  temperature  measurement,  there  is  no  need  to  view  the  entire  beam 
cross-section.  The  spatial  resolution  is  limited  only  by  the  volume 
viewed  by  the  optical  system.  However,  as  the  beam  spreads,  the  in- 
tensity of  the  luminescence  within  a volume  of  fixed  size  decreases. 
Hence,  for  temperature  measurements  the  main  effect  of  beam  spread  is 
to  reduce  the  signal -to-noise  ratio  of  the  measuring  system. 

The  shape  of  an  electron  beam  is  determined  by  the  interaction  of 
many  complex  phenomena,  particularly  at  high  densities.  For  very  low 
densities  and  short  path  lengths,  the  beam  electrons  can  be  assumed  to 
suffer  only  one  collision  as  they  traverse  the  field  of  interest. 
Theoretical  analyses  of  this  single  scattering  phenomena  can  be  conducted 
and  are  discussed  in  Ref.  1,  As  the  gas  density  is  increased,  both 
multiple  and  plural  scattering  occur  and  gas  focusing  effects  must  be 
considered. 

Low  energy  ions,  created  by  the  primary  excitation  process,  can  be 
trapped  within  the  potential  well  of  the  beam.  These  ions  can  reduce 
the  space  charge  repulsion  effects  within  the  beam  and  lead  to  phenomena 
generally  termed  gas  focusing.  Secondary  electrons  are  also  created  by 
the  excitation  process  and  their  presence  within  the  beum  and  the  sur- 
rounding region  may  have  a marked  effect  on  the  beam  shape.  Tlie 
secondary  electron  concentration  has  been  estimated  by  Harbour,  et  al, 
and  rather  large  electron  number  densities  were  predicted  at  certain 
operating  conditions.  Magnetic  pinch  effects  are  also  possible,  although 
these  are  expected  to  be  weak  for  electron  energies  typical  of  beams 
used  for  diagnostic  purposes. 

Gas  focusing  and  scattering  of  the  beam  electrons  within  the  beam 
generator  further  modify  the  beam  shape  in  the  test  gas.  Most  beam 
generators  employ  a small  exit  orifice  for  the  electrons  to  pass  from 
the  generator  into  the  test  gas.  A dynamic  pumping  system  is  employed 
to  maintain  the  generator  region  at  a sufficiently  low  pressure  (e.g., 
less  than  10”’  torr).  As  the  test  gas  pressure  is  raised,  scattering 
of  the  electrons  within  the  generator  occurs.  The  effects  of  this 
scattering  are  cancelled,  in  part,  by  adjusting  internal  deflection  and 
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focus  control  systems  to  maximize  the  beam  current  within  the  test  gas. 
At  high  densities,  maximization  of  the  beam  current  can  have  large 
influences  on  the  electron  trajectories  through  the  generator  and  in  the 
test  gas. 

Reliable  analytical  models  of  the  structure  of  the  beam-generated 
plasma  within  the  test  gas  are  not  available  because  of  the  extreme  com- 
plexity of  the  region.  Estimates  of  electron  trajectories  must  include 
the  effects  of  focusing  and  scattering  within  the  beam  generator,  in 
addition  to  the  mechanisms  active  in  the  high  density  test  region. 
Empirical  correlations  of  the  beam  shape  must  generally  be  used  to  de- 
scribe the  variations  of  beam  size  with  gas  density,  accelerating 
voltage,  and  path  length. 

Some  guidance  in  developing  empirical  correlations  can  be  obtained 
from  the  results  of  multiple  and  plural  scattering  theories  summarized 
in  Ref.  1.  The  scattering  can  be  characterized  in  terms  of  the  mean 
number  of  scatterings,  ft,  occurring  in  a length  L of  gas  with  a number 
density  of  ng,  a mass  m and  an  atomic  number  of  Z.  The  mean  number  of 
scatterings  is  given  by1 

ft  2.5  nZx/:’(Z  + l)ngI,hL7ra"v^  (34) 

where  ve  is  the  electron  velocity.  Data  summarized  in  Ref.  1 show  that 
the  angle  subtended  at  the  beam  origin,  ft,  by  one-half  of  the  character- 
istic lateral  spread  of  on  electron  beam  in  air  and  nitrogen  can  be 
represented  over  most  of  the  range  of  the  data  by 

ft  - C!iu  (35) 


where  c and  a arc  empirical  constants. 

Comae3'’  gives  correlations  of  the  beam  width  at  the  half- intensity 
point , Wj.  /2 , for  a 100  keV  beam  as 


’u: 


.044 


I/O 


(36) 


where  p is  in  torr,  I.  in  cm,  v«  in  kV,  T in  K,  and  W,  in  mm.  This 
can  be  put  in  the  form  of  Eq.  (35)  by  observing  that  the  angle  ft  is 
approximately  (Wv/t/2L).  Hence,  in  inks  units,  Eq.  (36)  becomes 

Wj.  /r?/L  = 6.6  x 10:'[pL/veT]1/1  (37  ) 

or,  combining  Eqs.  (3/,)  and  (37)  for  air, 

W,  ip. /Is  = .16  sfn  (38) 

Camac  includes  data  only  for  a beam  accelerating  potential  of  100  kV  and 
there  are  insufficient  data  to  judge  the  goodness  of  fit  of  Eq.  (36). 
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Wallace*  * gives  correlations  01'  W,  for  a 40  keV  beam  at  test  gas 
pressures  of  .85  and  1.6  Lorr.  Although  the  spatial  resolution  of 
these  measurements  was  relatively  poor  for  determining  beam  shape,  that 
portion  of  the  data  for  which  the  beam  size  was  larger  than  the  urea 
viewed  by  the  optical  system  agreed  reasonably  well  with  the  correlation 
of  Eq  . ( 16  ) . 

Detailed  experiments  were  conducted  in  the  static  chamber  to  deter- 
mine the  form  of  the  variation  of  Wv /n  over  a wide  range  of  pressure, 
beam  voltage,  and  beam  path  length.  The  accelerating  voltage  was  varied 
from  20  to  40  kV  with  air  pressures  from  0.2  to  30  torr.  The  width  of 
the  half-intensity  point  within  the  beam  was  obtained  by  scanning  across 
the  beam  with  the  spectrometer  set  to  view  the  head  of  the  N >+  first 
negative  (0,0)  bund  at  3914  K.  Hence,  the  half-width  data  reflect  only 
the  radiation  from  the  first  negative  system;  radiation  from  other 
systems,  particularly  the  N.->  second  positive  system,  was  not  included 
in  the  beam  shape  measurements.  This  will  affect  the  width  measurements, 
but  gives  data  most  meaningful  in  applications  of  an  electron  beam  to 
temperature  measurements  at  high  densities. 

The  experimental  data  are  presented  us  Wt  ,o/L  vs.  pi, /v^  in  Fig.  10. 
In  this  form,  Cumae's  equation  is 

Wws/L  .7!pft(torr)T,(in.)/vb(kV)]l/:'  (10) 

where  p should  be  interpreted  as  equivalent  pressure . The  fit  of 
Eq.  (39)  to  the  data  of  Fig.  10  is  relatively  poor,  due  primarily  to 
the  initial  beam  size. 

The  form  of  Kq.  (39)  can  be  justified  in  terms  of  multiple  scatter- 
ing theory  when  Vir;/h  is  taken  as  a measure  of  the  angle  of  the  beam 
electrons.  When  the  beam  spreading  is  small,  a better  measure  of  the 
angle  of  the  electrons  can  be  obtained  by  examining  the  change  in  the 
half -width  as  the  beam  tranverses  the  test  gas.  Accordingly,  the  half- 
width of  the  beam  at  the  exit  orifice  was  subtracted  from  the  measured 
half’-width.  The  data  are  presented  in  the  form  (Wt / .c-W-, ) /I , vs.  ph/vg 
in  Fig.  1; , where  W,,  is  the  minimum  beam  half  width. 

Equation  (39)  is  seen  to  give  a much  better  representation  of  the 
beam  shape  when  it  is  applied  in  terms  of  the  cnange  in  beam  width. 

The  scatter  in  Fig.  31  results  from  the  complex  interaction  of  the 
phenomena  causing  beam  spreading  which  cannot  be  represented  by  a simple 
scattering  theory.  In  addition,  the  effects  of  flow  on  the  beam  shape 
have  not  been  included  either  in  the  development  of  the  correlation  frern 
scattering  theory  or  in  the  experimental  measurements.  The  gas  flow 
will  remove  both  ions  and  secondary  electrons  from  the  beam  region, 
changing  the  focusing  effects  observed  in  the  static  case.  However, 
the  correlation  of  Eq.  (39)  is  adequate  for  estimating  the  beam  size  in 
air  and  nitrogen  test  gases  for*  both  stat  ic  r.nd  flov  cor.  lit. ions*. 
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V.  EFFECTS  ON  PREVIOUS  MEASUREMENTS 


The  studies  discussed  here  were  motivated,  in  part,  by  certain  dis- 
parities in  predicted  and  measured  nitric  oxide  number  densities  and 
vibrational  temperatures  in  previous  tests  conducted  in  the  AFFDL  Electro- 
gasdynamie  Facility  (EGF).  It  was  suspected  that  these  disparities  might 
result  from  inapplicability  of  the  usual  excitation-emission  model  used 
to  interpret  the  measurements6.  A summary  of  the  previous  data  obtained 
in  the  EGF  is  given  in  Table  V.  The  tests  conducted  here  with  nozzle  N2 
simulated  the  conditions  obtained  in  the  EGF  with  the  19  inch  nozzle. 
Evaluations  of  the  possible  effects  on  the  previous  data  resulting  from 
the  present  studies  are  given  below. 

A.  MOLECULAR  NITROGEN 

Hie  vibrational  temperature  and  number  density  of  molecular  nitrogen 
were  measured  in  the  EGF  with  both  the  19  inch  and  7 inch  nozzles;  the  N;> 
rotational  temperatures  were  measured  only  in  the  19  inch  nozzle.  As  seen 
in  Table  V,  Nj  rotational  temperature  data  were  obtained  at  equivalent 
pressures  below  approximately  0.25  torr.  At  these  equivalent  pressures  it  is 
well  known1  that  N?  rotational  tempera  tui*e  data  can  be  obtained  with  errors 
no  greater  than  5$. 

Tne  studies  conducted  here  demonstrate  that  N~  vibrational  tempera- 
tures can  be  measured  at  densities  at  least  an  order  of  magnitude  greater 
than  thess  employed  in  the  EGF.  Hence,  there  are  no  errors  introduced  in 
the  previous  vibrational  temperature  data  due  to  elevated  densities. 

The  accuracy  of  the  N2  number  densities  is  directly  related  to  the 
accuracy  of  the  quenching  constants  employed  to  reduce  the  data6;  however, 
since  the  maximum  equivalent  pressure  in  the  EGF  was  apiy_.  ximately  2 torr, 
relatively  large  errors  (30%)  in  the  quenching  constants  can  be  tolerated 
before  errors  greater  than  10$  in  ■ • .umber  doiisities  result.  Hence,  no 

additional  errors  in  N2  number  dens  result  from  the  present  studios. 

B.  MOLECULAR  a YGKij 

Data  for  molecular  oxygen  were  obtained  in  the  EGF  only  with  the 
19  inch  nozzle  where  the  equivalent  pressures  were  below  0.25  torr.  The 
present  studies  demonstrate  tnat  there  are  no  particular  anomolies  in  the 
0?+  second  negative  system  at  these  density  levels.  However,  under  the 
high  density  conditions  of  nozzle  N1  the  02  radiation  is  completely  over- 
lapped by  NO  > bands.  While  oxygen  data  for  the  7 inch  EGF  nozzle  .a 
not  yet  been  obtained,  it  is  likely  that  radiation  systems  other  than  u.  ■ 

Ci?+  second  negative  system  will  have  to  be  inves  tigated  for  use  in  measur- 
ing O2  '.'-ational  temperatures  in  high  density  air  flows. 
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C.  NITRIC  OXIDE 


Nitric  oxide  vibrational  temperature  and  number  density  data  were 
collected  in  the  EGF  with  both  the  19  inch  and  7 inch  nozzles  and  the 
data  were  interpreted  with  the  direct  excitation  model.  While  the 
results  of  the  present  studies  show  that  the  direct  excitation  model  is 
not  applicable  to  nitric  oxide,  the  use  of  the  scaled  band  intensity 
ratios  as  previously  discussed  eliminates  the  grass  error  in  band  inten- 
sity ratios  introduced  with  the  direct  model. 

As  can  be  seen  in  Fig.  26,  over  the  entire  density  range  investigated 
the  scaled  band  intensity  ratios  agree  as  well  with  the  direct  excitation 
mod.sl  as  they  do  with  the  vibrational  relaxation  model.  Neither  model 
adequately  describes  the  variation  of  band  intensity  ratio  with  vibrational 
temperature  and  excitation/de-excitation  mechanisms  not  presently  included 
in  these  models  may  be  important.  Because  of  the  uncertainties  in  the 
models,  there  is  at  present  no  justification  for  use  of  one  over  the  other 
so  long  as  the  measured  band  intensity  ratios  are  scaled  by  their  room 
temperature  values.  In  all  cases,  this  scaling  of  the  band  intensity 
ratios  was  employed  in  Ref.  6 so  that  further  corrections  to  the  data  are 
not  warranted. 


VI  . CONCLUSIONS 


A new,  systematic  study  of  electron  beam  induced  emission  has  bean 
conducted  at  density  levels  well  above  those  normally  used  in  applica- 
tion of  the  technique.  Emphasis  was  placed  on  the  fluorescence  induced 
in  air  flows,  since  there  is  much  interest  in  applying  the  electron  beam 
diagnostic  technique  to  flows  generated  by  arc-heated  wind  tunnels. 

'Ihe  results  of  the  wind  tunnel  studies  show  that  the  vibrational 
temperature  of  molecular  nitrogen  can  be  measured  with  good  accuracy 
at  densities  at  least  up  to  those  corresponding  to  an  equivalent  pres- 
sure^!' 20  €orr.  The  high  density  limit  for  application  of  the  diag- 
nostK  technique  to  N3  could  not  be  determined  from  the  present 
measurements.  Studies  should  be  performed  at  equivalent  pressures  at 
least.  ifco  100  torr  to  establish  this  high  density  limit. 

Comparisons  of  the  wind  tunnel  and  static  chamber  data  demonstrate 
that  the  static  chamber  is  unsuitable  for  examining  the  character  of 
electron  beam  induced  radiation  at  high  density.  Accumulation  of 
secondary  electrons,  low  energy  ions,  and  molecules  in  metastable  elec- 
tronic energy  states  contribute  to  the  excitation  of  the  ground  state 
molecules.  Such  low  energy  particles  are  not  present  in  meaningful  con- 
centrations in  flow  fields  with  high  velocities.  Similarly,  the  overlap 
of  the  N->'*  fiir.it  negative  system  by  the  N;.  second  positive  system  preva- 
lent in  the  static  chamber,  was  not  observed  in  the  wind  tunnel  sUidies. 


Since  it  is  known  that  the  second  positive  system  is  excited  primarily 
by  secondary  electrons,  this  is  further  indication  that  results  from  the 
static  chamber  bear  little  resemblance  to  those  obtained  at  the  same 
density  in  the  wind  tunnel  flow  field.  This  is  not  meant  to  be  a con- 
demnation of  all  previous  static  chamber  results.  The  degree  of  excita- 
tion by  secondary  particles  will  undoubtedly  be  a complex  function  of 
the  mechanical  and  electrical  arrangement  of  the  static  test  chamber. 
However,  such  experiments  must  be  conducted  with  great  care  and  the 
data  must  be  carefully  analyzed  to  remove  extraneous  effects  which  may 
not  be  present  in  high  speed  flow  fields. 

Comparisons  of  theoretical  and  measured  band  intensities  for  the 
NO  y system  show  that  the  direct  excitation  model  does  not  lead  to  the 
proper  relative  intensities  within  a given  band  sequence.  A new  theo- 
retical model  of  the  excitation-emission  process  has  been  formulated  to 
allow  for  transitions  between  vibrational  energy  levels  in  the  excited 
electronic  energy  state  of  tfte  NO  y system.  A relative  v'  - 1 -<■  v'  = 0 
transition  probability  was  determined  from  static  tests  and  applied  to 
the  wind  tunnel  data.  Comparisons  of  the  theoretical  and  experimental 
data  show  that  the  model  properly  accounts  for  the  general  trend  of  band 
intensity  ratio  with  vibrational  temperature  but  that  more  detailed 
studies  are  required  to  substantiate  the  validity  of  the  model. 

Wind  tunnel  and  static  chamber  experiments  were  conducted  with  the 
O',  second  negative  system.  At  equivalent  pressures  above  4.0  torr,  the 
observed  0a+  radiation  was  severely  overlapped  by  the  NO  y system. 
Sensitivity  studies  conducted  in  the  static  chamber  show  that  the 
expected  0a+  radiation  intensity  was  well  above  the  minimum  detectable 
level.  The  results  indicate  that  the  0->+  radiation  is  more  severely 
quenched  taan  is  that  in  the  NO  y system.  Hie  spectral  proximity  of 
the  0a+  and  NO  y bands  coupled  with  the  decreased  relative  intensities 
of  the  0a+  bands  at  high  density  cause  the  NO  y system  to  completely 
overlap  the  O,.'1  system  at  the  spectrometer  resolution  required  for  Ihese 
experiments . 

Beam  shape  experiments  were  conducted  in  a static  chamber  to 
establish  an  empirical  correlation  of  beam  half-width  with  density  and 
acceleration  potential.  An  analysis  based  on  the  angle  of  divergence 
of  the  beam  leads  to  an  acceptable  form  of  the  correlation. 
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